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Context. Complex turbulent motions are ubiquitously observed in many astrophysical systems. The origin of this turbulence, 
however, is still poorly understood. 

Aims. When cosmic structures form, they grow in mass via accretion from their surrounding environment. We propose that 
this accretion is able to drive internal turbulent motions in a wide range of astrophysical objects and study this process in the 
case of galaxies, molecular clouds and protoplanetary disks. 

Methods. We use a combination of numerical simulations and analytical arguments to predict the level of turbulence as a function 
of the accretion rate, the dissipation scale, and the density contrast, and compare with observational data. 
Results. We find that in Milky Way type galaxies the observed level of turbulence in the interstellar medium can be explained 
by accretion, provided that the galaxies gain mass at a rate comparable to the rate at which they form stars. This process is 
particularly relevant in the extended outer disks beyond the star-forming radius. In order to drive turbulence in dwarf galaxies, 
the accretion rate needs to exceed the star formation rate by a large factor and we expect other sources to dominate. We also 
calculate the rate at which molecular clouds grow in mass when they build up from the atomic component of the galactic gas 
and find that their internal turbulence is likely to be driven by accretion as well. It is the very process of cloud formation that 
excites turbulent motions on small scales by establishing the turbulent cascade. In the case of T Tauri disks, we show that 
accretion can drive subsonic turbulence at the observed level if the rate at which gas falls onto the disk is comparable to the 
rate at which disk material accretes onto the central star. This also explains the observed relation of accretion rate and stellar 
mass, M oc M*' 8 . The efficiency required to convert infall motion into turbulence is of the order of a few percent in all three 
cases. 

Conclusions. We conclude that accretion-driven turbulence is a universal concept with far-reaching implications for a wide range 
of astrophysical objects. 
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1. Introduction 

Astrophysical fluids on virtually all scales are character- 
ized by highly complex turbulent motions. This ranges 
from the gas between galaxies to the interstellar medium 
(ISM) within them, as well as from individual star- 
forming molecular clouds down to the protostellar ac- 
cretions disks that naturally accompany stellar birth, 
and has far reaching consequences for cosmic struc- 
ture formation. For example, it is the complex inter- 
play between supersonic turbulence in the ISM and self- 
gravity in concert with magnetic fields, radiation, and 
thermal pressure tha t determines when and where stars 
form in the Galaxy dLow fc Klessen 2004 Larson 2005 



Similar is true for protostellar accretion disks, where tur- 
bulent motions cause angular momentum redistribution 
and thus determine the rate at which material accretes 
onto the central star and the likelihood to build up plan- 
ets and planetary systems. 

Yet, despite its ubiquity and importance, very little is 
known about the origin of astrophysical turbulence. The 
number of possible sources is large and varies strongly de- 
pending on the physical scale under consideration. For a 
discussion of possible so urces of ISM turbulence, s ee e.g. 
Low fc Klessenl (|2004T ) or lElmeereen fc Scald (|2004l l. Here 



Ballesteros-Paredes et all 120071 : iMckee fc Ostrikerl (20071 ). 



we attempt to argue that it is the accretion process, that 
inevitably goes along with any astrophysical structure for- 
mation, let it be the birth of galaxies or stars, that drives 
the observed turbulent motions. We propose that this pro- 
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cess is universal and makes significant contr ibutions to 



the t urbulent energy on all scales (see also iField et al 



2008). We ask: Does the accretion flow onto galaxies, onto 
dense clouds in the ISM within these galaxies, and finally 
onto the protostellar accretion disks that accompany stel- 
lar birth within these clouds provide enough energy to 
account for the observed internal motions? What is the 
expected efficiency for the conversion of kinetic energy as- 
sociated with the infalling material into kinetic energy as- 
sociated with internal turbulence? Our analysis leads us 
to believe that accretion is indeed an important driver of 
turbulence on all scales observed. 

We structure our discussion as follows: We introduce 
the concept of accretion-driven turbulence in Section [2 
We estimate the energy input associated with accretion 
and compare it to the energy needed to compensate for 
the decay of turbulent motions assuming overall steady 
state. Under typical conditions the energy gain from accre- 
tion exceeds the energy loss by the decay of turbulence by 
far. However, we do not know the efficiency with which in- 
fall motions are converted into random turbulent motions. 
To get a handle on this quantity, we resort to numerical 
simulations of convergent astrophysical flows for guidance 
and propose a theoretical explanation of the trends in- 
ferred from the simulations. In Section [3] wc apply our 
method to galactic scales and propose that the turbulent 
velocity dispersion measured in the disk of the Milky Way 
and other galaxies are caused by accretion streams that 
originate in or pass through the halo. We then turn to 
the scales of individual interstellar gas clouds and argue 
in Section [4] that it is the process of cloud formation that 
drives their internal turbulent motions. Our third applica- 
tions lies on even smaller scales. In Section [5] wc speculate 
about the origin of turbulence in accretion disks. We focus 
our discussion on protostellar accretion disks during the 
late stages of the evolution (class 2 and 3 phases), but we 
note that similar arguments may apply to the accretion 
disks around black holes in active galactic nuclei. Finally, 
we conclude in Section [5] 

2. Basic Concept 

2.1. Energy Balance 

Several numerical studies ( Low et al.l 1998 ; Stone et al. 
19981: IPadoan k Nordlundl [19991 iLowl Q999I: lElmeereen 



20001) have demonstrated that supersonic turbulence de- 
cays on a timescale that is equivalent to the turbulent 
crossing time, 



a 



(1) 



where is the driving scale and a is the 3-dimensional 
velocity dispersion. This holds regardless whether the gas 
is magnetized or not and also extends into the subsonic 
regime. The exact value of Ld is not well constrained by 
the observational data and needs to be chosen with care 
for each system under consideration. We note, however, 



that it is a universal feature of the objects we study here 
that the bulk of the kinetic energy is carried by the largest 
spatial modes, consistent with turbulence being driven 
from the outside (see, e .g. lOssenkopf k Lowj|2002l : iBruntl 
2003 : Brunt et al. 20091 for nearby molecular clouds). 

The total loss of turbulent kinetic energy, E = 
1/2M<7 2 , to a system with total mass M through tur- 
bulent decay sums up to 



E, 



decay 



E 

Td 



1 Ma 3 
2~Td~ 



(2) 



When the system accumulates mass at a rate M the 
associated kinetic energy is 



1 . „ 
E la = -M in v in 



(3) 



where «i n is the infall velocity. 

We introduce an efficiency factor 



E 



decay 



(4) 



which represents the fraction, e, of the available accretion 
energy required to sustain the observed turbulent veloci- 
ties. 

For the hypothesis of accretion-driven turbulence to 
work, clearly E- m > -Edccay is required. This is usually 
true as we discuss in the Sections below. Wc note, how- 
ever, that the fraction of the infall energy that actually 
is converted into random turbulent motions is very diffi- 
cult to estimate. Clearly some fraction of the accretion 
energy turns into heat and is radiated away. In addition, 
if the system is highly inhomogeneous with most of the 
mass residing in high-densities clumps with low volume 
filling factor, most of the incoming flux will feed the ten- 
uous intcrclump medium rather than the dense clumps, 
and again, not contribute directly to driving their inter- 
nal turbulence. This is taken into account in the efficiency 
factor. Numerical experiments indicate that e depends on 
the density contrast between the infalling gas and the ma- 
terial in the system under consideration (see Section |2~2"|) . 
For molecular clouds forming in convergent flows e is of 
order of 0.01 to 0.1. If these values are representative for 
other systems, then in general E ln needs to be 10 to 100 
times larger than -Edecay 

2.2. Estimate of Efficiency 

To estimate the efficiency at which accretion energy is 
converted into turbulent energy we resort to numerical 
simulations of converging f lows ( e .g. lAudit k Hennebelle 
20051: iHeitsch et al.1 [2005 1 l2006at iFoiini k Warder! 120061: 



200 



Vazauez-Semadeni et al.l 20061 2007 : Hennebelle et al 



Audit fc Hennebelle! l2010l ; iBaneriee et al.l l2009n T 



These simulations consider two colliding flows of diffuse 
gas which produce strong density fluctuations of cold gas. 
The incoming velocity is initially supersonic with respect 
to the cold and dense gas which forms under the influence 
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of cooling and ram pressure. It is generally found, that the 
resulting turbulence in this component is comparable to 
the observed values in Galactic molecular clouds. 

The simulations rep orted here are v ery s imilar 
to those presented by Hennebelle et al. ( 20081) and 
Audit fc Hennebelle! (|2010h . Thev have been performed 



with the adaptive m esh-refinement magnetohydrodynam 



ics code RAMSES (jTevssierl 120021 iFromang etaH I2OO6I) 



and include magnetic fields (initially uniform and equal 
to 5 /iG) and self-gravity. They are either isothermal at 
T = 50 K or start with a warm neutral medium at 
T = 8000 K and self-consistently treat cooling processes 
assuming a standard 2-phase ISM cooling function. The 
gas is injected from the boundary with a density equal 
to 1cm -3 and a mean velocity of either 15 or 20kms~ 1 
on top of which fluctuations with an amplitude of 50% 
have been superimposed. In order to quantify the impact 
of the numerical resolution (a crucial issue) as well as the 
influence of the thermal structure of the flow, we present 
the results of four calculations. First, two lower resolution 
calculations with an incoming velocity of 15kms _1 , one 
isothermal and one in which cooling is treated. Both have 
an initial grid of 256 3 computing cells and two further 
AMR levels are used when the density reaches a threshold 
of 80 and 160 cm -3 . Second, we also present two simula- 
tions with a higher incoming velocity of 20 km s" 1 . One 
has a high resolution and starts with 512 3 computing cells 
and four further AMR levels arc used when density reaches 
respectively 50, 100, 400 and 1600 cm -3 . In this calcula- 
tion the number of cells is about 5 x 10 8 . The other has the 
same resolution as the two lower resolution simulations. 
Figure [T] shows the column density in the computational 
box for the high resolution run. 
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Fig. 1. Column density at t = 18.75 Myr in the high res- 
olution colliding flow calculation. 



Fig. 2. Mass, velocity dispersion and efficiency of the en- 
ergy injection as a function of gas density in four collid- 
ing flow calculations. Solid, dotted and dashed lines show 
the cases with standard ISM cooling. The solid one corre- 
sponds to the highest numerical resolution and an incom- 
ing velocity of 20kms~ 1 , the dash-dotted line is identical 
except that it has a lower resolution, and the dotted line 
is for a lower resolution simulation and an incming veloc- 
ity of 15kms _1 . The dashed-dotted line corresponds to 
the purely isothermal calculation (with gas temperature 
of about 50 K). It has the same resolution than the lower 
resolution simulations with cooling. 

Figure [5] shows the efficiency e as defined by Eq. 
as a function of the gas density. That is, we select all 
cells in the computing domain with densities above a 
certain threshold value pt and then compute the quan- 
tity J pv 3 dV = M t (7 3 , where p and v are the density and 
velocity of the cells, and M t and er t are total mass and ve- 
locity dispersion of the gas above p t ■ Finally we divide by 
Lbox and by Ei n as defined by Eq. ([3]). Various trends can 
be inferred from Fig. [2] First, for all the simulations we 
find that the efficiency decreases with the gas density as 
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roughly 1/p. The two low resolution simulations with cool- 
ing but different incoming velocities are very close to each 
other. The high resolution simulation exhibits a slightly 
higher efficiency. This is expected as the numerical dis- 
sipation is lower. However, it is larger by only a small 
factor of 1.5 to 2, suggesting that our result is reason- 
ably well converged. Finally we see that the efficiency is 
larger in the isothermal case by a factor of ~ 3 compared 
to the low resolution runs with cooling. This indicates 
that the thermal structure of t he flow has a significant im - 
pact on its dynamics (see also Audit fe Hennebelle 201dt ). 
Because the isothermal run and the simulations with a 2- 
phase medium cooling function are significantly different, 
the discrepancy gives us an estimate of the uncertainty 
of the efficiency. Interestingly, Fig. [5] shows two different 
regimes. At low densities the total mass above the thresh- 
old value decreases slowly with density, while the veloc- 
ity dispersion decreases more steeply. At large densities, 
however, the mass decreases rapidly with density while 
the velocity dispersion is nearly constant. This behavior 
is discussed further in Appendix iBl 

We conclude that for astrophysical systems, such as 
molecular clouds, of mean density p accreting gas at den- 
sity ~ pi n , we typically have 



Pin/P- 



(5) 



This relation is expected to be valid within a factor of a 
few. 



2.3. Possible Explanation for the e - p 1 Relation 

As the e - p _1 relation appears to be both important and 
interesting, we propose a possible theoretical explanation. 
Consider a turbulent flow with wave numbers ranging from 
fcmin to fc max . For incompressible fluids, the power spec- 
trum of the velocity field based o n dimensional argu ments 
is expected to be E(k) oc k~ b / 3 (jKolmogorovl 1 1 94 11 ) . The 
kinetic energy carried by wave numbers k and larger is 
given by the integral j^ max E(k)dk, with the total kinetic 



energy being J fc max E(k)dk. For compressible media, the 
scaling relation is more complica ted as the density de - 



pendency needs to be considered (jvon Weizsackerlll951l ). 
Based on the argument that (pv 3 /l)/e e is dimensionless 



dFerrini et alj 


19831: Fleck 1983. 1996: iKritsuk et al.1 


2007: 


Schmidt et al. 


2008) that the relation E(k) oc k~ 5/3 still 



holds, provided that E{k) is the power spectrum of p x l 3 v 
instead of v. 

Density fluctuations in a turbulent flow follow a 
roughly log-normal behavior. When identifying clumps 
and cores, e.g. defined as connected groups of cells/pixels 
above some thresholds, it has been found that their mass 
spectrum often follows a po wer law dN / dM oc M a with a 



or 



— 1. 7 (see, e.g. iHeithausen et al.lll998| for the 



slope of a ps 
obser vations 
2006L 

tions, and iHennebelle fc Chabrie 



~~ — II 

Klessen 2001, Ballest cros-Paredes et al 
Hennebelle fc Auditl 120071 for n umerical simula- 
ierl 12008 for analytical ar- 



guments). In a convergent flow the biggest clumps cannot 
be much larger than l t = (pt/po) 1 x Lq where po and Lq 
are the typical density and scales of the large scale flow. 
Thus, we expect that the largest scale at which clumps 
denser than p t exist, is it oc p^ 1 . 

The quantity f?J l 'a 2 integrated over the cells denser 
than pt, is thus expected to be of the order of 
j-^max fiQ^flfc^ where E is the power spectrum of p 1 / 3 v 



and kt — 2n/lt- Thus, it is found that (p 2 ^ 3 of ) 
f t /3 oc p" 2/3 . This leads to (p 2 t /3 a 2 ) 3 / 2 « (p t a 3 ) oc p t " ] 
and after multiplication by the volume of the cloud 
M t a 3 cxpt 1 . 



!.-2/3 



So far, we have simply shown that the quantity M t a 3 
obtained by integration over scales smaller than l t is pro- 
portional to it but it could be the case that the dense 
parts of the gas, i.e. regions denser than p t , have a negli- 
gible contribution to this integral, in particular because of 
their low filling factor. However two arguments are in dis- 
favor of this statement. First, Hennebelle & Audit (2007) 
have calculated the power spectrum of the kinetic energy 
of the flow, pv 2 , while clipping dense structures above var- 
ious threshold values. As can be seen in their Fig. (14), 
the energy contained in large-scale motions is unchanged 
when varying the threshold whereas the energy contained 
on the small scales (below about one hundredth of the 
computing box length) decreases with increasing density 
threshold. It is dominated by high-density structures. The 
second argument is also inferred from numerical simula- 
tions. The power spectrum of v has been calculated in sev- 



eral studies (e.g. iKlessen et al. I l200d; iHeitsch et al. I [200! 



Kritsuk et all 120071: iFederrath et al.1 120091 ) and typically 
it has been found to be P(v) oc fc -19 . Thus, considering- 
only gas at densities close to the mean value of the sys- 



tem, p, we infer that J fe max p 2 l 3 v 2 dk oc k t 9 oc I 



°- 9 while 



Jfct 



p 2 l 3 v 2 dk 



-2/3 



2/3 



The implication is that 



as It decreases, the contribution of the energy contained 
in scales smaller than l t due to the diffuse gas becomes 
smaller and smaller with respect to the energy contained 
in the dense gas at these scales. 

The relation, M t of oc p^ 1 , is therefore broadly con- 
sistent with the trend we measure. The coefficient seems 
more difficult to predict and as our numerical simulations 
suggest, it may vary from one flow to another. This re- 
quires further investigation. We note in this context that 
the situation is strongly reminiscent of purely incompress- 
ible turbulence where dissipation occurs in a subset of 
space, in filaments with small filling factor but high vor- 
ticity. It can be described by a multi-f ractal statistical 
approach (jFrisch et al.lll978t lFrischlll995l ) to take into ac- 
count their intermittent nature. For a model based on en- 



ergy dissipation in shock-generated sheet s, see iBoldvrev 



1994) 



(2002) extending the theory developed bv lShe fc Levecme 
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3. Turbulence in Galactic Disks 

3.1. General Considerations 

In our first application we investigate the question as to 
whether accretion from an external gas reservoir could 
drive the velocity dispersion observed in spiral galaxies. 

The Milky Way, as a typical galaxy, forms new 
stars at a rate of M$f ~ 2 — 4 M w yr~ 1 . Its gas mass 
out to 25 kpc i s - 9 x 10 9 M o (|Naab fc Ostrikerl 120061; 
Xue et aT]|2008 ). Assuming a constant star formation rate, 
the remaining gas should be converted into stars within 
about 2 — 4 Gyr. Similar gas depletion timescales of order 
of a few billion years are rep orted for many nearby spiral 



galaxies ( Bigiel et al. 20081) . This is much shorter than 



the ages of these galaxies which is ~ 10 10 yr. If we discard 
the possibility that we observe them right at the verge 
of running out of gas, and instead assume they evolve in 
quasi steady state, then these galaxies need to be supplied 
with fresh gas at a rate roughly equal to the star formation 
rate. 

The requirement of a steady accretion flow onto typical 
disk galaxies is a natural outcome of cosmological struc- 
ture formation ca lculations if baryon ic ph ysics is consid- 
ered c onsistently. iDekel et al.1 (|2009h and ICeverino et al 



( 20091 ). for example, argue that massive galaxies are con 
tinuously fed by steady, narrow, cold gas streams that pen- 
etrate through the accretion shock associated with the 
dark matter halo down to the central galaxy. Roughly 
three quarters of all galaxies forming stars at a given rat e 
are fed by smooth streams (see also Agertz et al. 20091 ). 
On large scales, also the fact that the observed amount 
of atomic gas in the universe appears to be roughly con- 
stant since a redshift of z sa 3 although the stellar content 
continues to increase, suggests that HI is continuously re - 
plenished ([Hopkins et al.ll2008t Iprochaska fc Wolfell2009l ). 
For our Galaxy, further evidence for a ongoing inflow of 
low-metallicity material comes fro m the presence of deu- 
terium at the solar neigh borhood (|Linskyl 2003 ) as well 
as in the Galactic Center ( Lubowich et al.ll2000 ). As deu- 
terium is destroyed in stars and as there is no other known 
source of deuterium in the Milky Way, it must be of cosmo- 
logical and extragalac tic origin (jOstriker fc Tinslevlll975 : 
Chiappini et al.ll2002l ). 

It is attractive to speculate that the population of 
high-velocity clouds (HVC) observed around the Milky 
Way is the visible signpost for high-density peaks in 
this accretion flow. Indeed the inferred HVC infall rates 
of 0.5 - SM^vr- 1 (jWakker et alJll99flt iBlitz et al.lll999l 
Braun fc Thilkerl 120041 |Putmanll20o'i r are in good agree- 
ment with the Galactic star formation rate or with 



chem ical enrichment models (see, e.g. ICasuso fc Beckman 
20041 and references therein). An important question 
in this context is where and in what form the gas 
reaches the Galaxy. Recen t numerical simulations indicate 
( Heitsch fc Putmanll2009l ) that small clouds (with masses 
less then a few 10 4 M Q ) most likely will dissolve, heat up 
and merge with the hot halo gas, while larger complexes 



will be able to deliver cold atomic gas even to the inner 
disk. We explore the idea of continuous gas accretion onto 
galaxies and argue that this process is a key mechanism 
for driving interstellar turbulence. 

One of the remarkable features of spiral galaxies is 
the nearly constant velocity dispersion er, e.g. as mea- 
sured in HI emission lines, rega r dless of galaxy mass and 
type (iDickev fc Lockman 1990t van Zee fc Bryant 19991 
Tamburro et all 120091) . The inferred values of a typi- 
cally fall in a range between 10 km s" 1 and 20kms~ 1 
(|Bigiel et al.ll2008t IWalter et al.ll2008l) and extend well be- 
yond the optical radius of the galaxy with only moder- 
ate fall-off as one goes outwards. It is interesting in this 
context that the transition from the star-forming parts 
of the galaxy to the non-star-forming outer disk seems 
not to cause significan t changes in the velocity dispersion 
( Tamburro et al. 2009J). This apparent independence from 
stellar sources sets severe constraints on the physical pro- 
cesses that can drive the observed level of turbulence. 

Sever al possibilities have been discussed in the lit 



erature dLow fc Klessen 



2004; lElmegreen fc Sca lo 2004; 



Scalo fc Elmegreenll2004l ). Large-scale gravitational insta- 
bilities in the disk, i.e. spiral de nsity waves, can potentially 



provide sufficient energy (e.g. iLi et al.ll20"05h . However, 



the efficiency of this process and the details of the cou- 
pling mechanism arc not well understood. Similar holds 
for the magneto-rotational instability (MRI) which has 
been identified as a main source of turbulenc e in proto- 
stellar accretion disks ( Balbus fc Hawlevlll998l ). Although 
we have amp le evidence of the presence of large-s cale mag- 
netic fields (|Heiles fc Trolandl l2005t iBeckl 120071) there is 
some debate whether the MRI can provide eno ugh energy 



to explain tne observed levels ot turbulence imecK . et ai. 
19961 ISellwood fc Balbuslll999l: iDziourkevitch et al.ll2004l: 
Piontek fc Ostrikerl 20071 ). For the star-forming parts of 



spiral galaxies, clearly stellar feedback in form of expand- 
ing HII bubbles , winds, or supe r nova explosions plays an 
important role. iLow fc Klessen! ( 2004 ) show that the en- 
ergy and momentum input from supernovae is a viable 
driving mechanism for interstellar turbulence. However, 
this approach clearly fails in the extended outer HI disks 
observed around most spiral galaxies and it also fails in 
low-surface brightness galaxies. Here accre tion driven tur 
bulence seems a viable option (see, e.g. 
20071) . 



Santillan et al 



3.2. Energy Input Rate 

In order to calculate the energy input rate from the ac- 
cretion of cold gas we need to know the velocity v ln with 
which this gas falls onto the disk of the galaxy and the effi- 
ciency e with which the kinetic energy of the infalling gas 
is converted into ISM turbulence. As the cold accretion 
flow originates from the outer reaches of the halo and be- 
yond and because it lies in the nature of these cold streams 
that gas comes in almost in free fall, i>; n can in principle 
be as high as the escape velocity u eS c of the halo. For the 
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Milky Way in the solar neighborhood v, 



550 km; 



( Fich fc Tremaindll99lt ISmith et al.ll2007l) . However, nu- 
merical experiments indicate that the inflow velocity of 
cold s treams is of order of the virialization velocity of the 
halo (jDekel et al.ll2009( ) which typically is ~ 200 km s" 1 . 
The actual impact velocity with which this gas interacts 
with disk material will also depend on the sense of rota- 
tion. Streams which come in co-rotating with the disk will 
have smaller impact velocities than material that comes in 
counter-rotating. To relate to quantities that are easily ob- 
servable and to within the limits of our approximations we 
adopt Wi n = v rot as our fiducial value, but note that con- 
siderable deviations are possible. We also note that even 
gas that shocks at the virial radius and thus heats up to 
10 5 — 10 6 K, may cool down again and some fraction of it 
may be available for disk accretion. This gas can condense 
into higher- density clu mps that s ink towards the center 
and replenish the disk ( Peek! 2009h . Again, V[ n ps v mt is a 
reasonable estimate. 

We can now calculate the energy input rate associated 
with this accretion flow as 



E„ 



= +1.3 x 10 40 ergs _1 • 

( Mn 



1 M© yr- 



200 kms" 



(6) 



By the same token, the energy loss through the decay of 
turbulence is 



E. 



decay 



E 



1 Ma 3 



3.2 x 10 39 ergs _1 
M \ ( a 



10 9 M, 



lOkms- 



100 pc 



(7) 



Given an efficiency e as defined in Eq.|4]the mass accretion 
rate required to compensate for this energy loss is 



1 2E, 



decay 



0.25 M yr" 1 



1 Ma 3 
1 



M 



10 9 Mr. 



10 km; 



La 



100 pc 



200 kms 



(8) 



3.3. The Milky Way 



Current mass models of the Milky Way (jXue et al.ll2008l) 
indicate a total mass including dark matter of about 
1 x 10 12 M Q out to the virial radius at ~ 250 kpc. The 
resulting rotation curve is 220 km s -1 at the solar radius 
Rq ps 8.5 kpc and it declines to values slightly below 
200 km s -1 at a radius of 60 kpc. The total mass in the 
disk in stars and cold gas is estimated to be ~ 6 x 10 10 M Q . 
Assuming a global baryon fraction of 17% this corresponds 
to 40% of all the baryonic mass within the virial radius 



Table 1. Properties of gas components of the Milky Way. 

Component molecular gas atomic gas e 
M (10 9 M o ) a 
L d (pc) b 
er (kms^ 1 ) 1 
E kin (10 55 ergs- 1 ) d 



2 


6 


150 


1000 / 


5 


12 


0.5 


8.6 



a T otal mass of the component. Values from [Ferrierel ((2001) 
and iKalberlal l|2003l ). 

b We take twice the observed disk scale height as the true thick- 
ness to be considered. 

c The parameter a is the 3-dimensional velocity dispersion. We 
take the 1-dimensional velocity dispersion of the molecular gas 
to be 2.9kms~ 1 and of the HI gas to be 6.9kms _1 . 
d Total kinetic energy of the component, -Ekin = 1/2 Ma 2 . 
e The atomic component in principle can be separated into a 
cold (T ps few x 10 2 K) and a hot (T ps few x 10 2 K) compo- 
nent. Because they have similar overall distribution we consider 
them together and take - whenever possible - mean values. 
' The scale height of HI ranges from ~ 230 pc within 4 kpc 
up to values of ~ 3 kpc at the outer Galactic boundaries. The 
HI disk therefore is strongly flared. We adopt some reasonable 
mean value, but note that this introduces additional uncer- 
tainty. 



and implies that roughly the same amount of baryons is 
in an extended halo in form of hot and tenuous gas. The 
gaseous disk of the Milky Way can b e decomposed i nto a 
number of different phases. We follow Ferrierel ( 2001 ) and 
consider molecular gas (as traced, e.g. via its CO emis- 
sion) as well as atomic hydrogen gas (as observed, e.g. by 
its 21cm emission). The adopted values are summarized 
in Table Q] We neglect the hot ionized medium in our 
analysis, as Galactic HII regions are produced and heated 
predominantly by the UV radiation from massive stars 
and thus should not be included here. We note that 95% 
of the turbulent kinetic energy is carried by the atomic 
component. 

If we use the numbers from Table [TJ assume M- m = 
Msf w 3M Q yr _1 , and adopt the fiducial value v- m — 
220 kms" 1 , then Eq.'s © and © yield 



E in ps +4.6 x lO^ergs" 1 , 
Sdccay PS -3.9 x 10 39 ergs" 1 , 

requiring an efficiency of only 

£= |^docay|/^in W0.08. 



(9) 
(10) 



(11) 



In the light of Eq. ([51), thi s is a r easonable number. If we 
follow Heitsch fc Putman (2009) and adopt densities in 
the range n = 0.01 to 0.1 enr 3 for the accreting gas clouds 
and assume a mean ISM de nsity of 1 cm- 3 in the solar 
neighborhood (|Ferrierel200ll) as well as a drop to 0.1cm 3 
out at a distance of 25 kpc, we expect the efficiency to be 
of order of 10%. 
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3.4. Spiral Galaxies 



The HI Nearby Galaxy Survey, THINGS, (jWalter et al 
2008) opens up the possibility to perform the above anal- 
ysis for the extended HI disks of other spiral galax- 
ies as well. We obt ained the dataset discussed by 
Tamburro etafl (|2009t) . which allows us to analyze HI 
column density Shi and vertical velocity dispersion ctid 
as function of radius for 11 nearby galaxies. Our sample 
includes 8 Milky Way type spirals as well as 3 gas-rich 
dwarfs. Each galaxy map contains between 140.000 and 
720.000 data points. The main parameter s are provided 
in Table [21 Please consult Tamburro et al. ( 20091 ) for fur- 
ther details on the original data set. 

We integrate over the entire map to get the to- 
tal HI mass, and convert ctid to the 3-dimcnsional ve- 
locity dispersion cr assuming isotropy. We read off the 
measured star for mation rate, Msf, from Table 1 of 
Walter et al. I (l2008h. and consider the rotation curves from 



de Blok et al.l ( 20081 ) , where we adopt the peak value v rot 
for our analysis. An estimate of the tur bulent length scale 
id is more difficult to obtain. We follow lLerov et al.1 (|2008j) 
and derive an estimate of the thickness H of the HI layer 
by assuming hydrostatic equilibrium at every pixel. In this 
case, 

H « a 2 1D /(2irGZ m ) , (12) 

with cr id = 3~ 1//2 cr and G being the 1-dimensional veloc- 
ity dispersion and the gravitational constant, respectively. 
Note, that we neglect the contribution from molecular gas 
as well as from stars in the determination of H which could 
be significant especially in the inner regions of the disk. In 
principle we would need to calculate H ps a- 2 D /(27rGI]tot) 
with S to t = Shi + Sh 2 + S* being the combined surface 
density of gas and stars. Our estimate of the scale height H 
in the inner parts of the galaxy therefore is an upper limit. 
An alternative estimate for the disk thickness is based on 
the total enclosed mass at any given radius R using the 
local circular velocity v TOt (see Appendix |A"]) , 



H « i?cr 1D /v rot . 



(13) 



Again, we follow lLerov et al.l (|2008l ) and approximate the 
rotation curve with the fit formula, 



Vrot (R) = Wflat 



1 — cxp 



R 

-Rflal 



(14) 



using the values ffl a t and i?fl a t for the flat parts of the ro- 
tation curve from their Table 4. We point out, that except 
maybe for the inner disk the mass distribution is domi- 
nated by the dark matter content of the galaxy. 

Once we have calculated H at each location in the 
map, we set the local turbulent scale length to id = 2 H. 
For most galaxies in the sample, both of the above es- 
timates lie within a factor of two or less of each other, 
with the potential method usually giving somewhat lower 
numbers for the large spirals and higher values for the 
dwarf galaxies. The exceptions are NGC 5193, NGC 3351, 
NGC 5055, and NGC 4736 which have a strong molecular 



component in the center ( Lerov et al.ll2008( ) and where our 
id estimate based on vertical hydrostatic balance using 
Shi consequently is too large. NGC 5055 and NGC 4736, 
furthermore, are characterized by extended HI streamers 
at ~ i?2 5 which results in a locally enhanced velocity dis- 
persion ( Walter et al.l[2008h . again leading to inflated id 
values from hydrostatic balance. For comparison with the 
Milky Way, mean velocity dispersion (cr), mean HI sur- 
face density (Shi), and mean turbulent scale (id) based 
on both methods are obtained for each galaxy as surface 
density weighted average over all pixels and are provided 
in Table Has well. 

The radial variation of cr, Shi, and id is shown in 
Figure [3l In addition, this figure also indicates the num- 
ber of pixels contributing to radial annuli of width AR = 

0. 5 kpc for each of the THINGS maps. This provides some 
estimate of the statistical significance of the data at differ- 
ent radii. We notice that outside of R25, which is a proxy 
for the optical radius of the stellar disk, the values of cr 
typically drop below ~ 20 kms" 1 and can get as low as 
~ 5 kms™ 1 in the outer disk of the dwarf galaxies in the 
sample, with statistical uncertainties of 2—3 kms -1 . Inside 
of R25 the velocity dispersion can be significantly higher. 
Here stellar source s provide additional e nergy for driving 
turbulent motions ( Low fc Klessenl 20041 ) . The HI surface 
density also drops outside of i? 2 5- We note, however, that 
for some galaxies with very extended disks, Shi can re- 
main as high as lM Q pc -1 out to 4i? 2 5- We also point 
out, that some galaxies reveal a noticeable HI depletion 
in t heir central regions wh ere the gas is mostly molecu- 
lar dLerov et al.ll2008ll2009l ). Many galaxies in our sample 
show significant flaring in the outer disk. Some galaxies 
also show large id values in the very inner parts. Recall 
that this is in part an artifact of our analysis, because 
we neglect the contribution to the graviational potential 
from the stellar disk and from the molecular gas when 
computing the scale height from vertical hydrostatic bal- 
ance. Inflated central id therefore correlate well with the 
inner depletion of Shi- 

With this information, we can now compute the local 
turbulent kinetic energy as well as the local kinetic en- 
ergy decay rate based on the two estimates of the disk 
scale height. We display the result in the left two columns 
of Figure |3J For simplicity, we only show the decay rate 
based on the assumption of vertical hydrostatic equlib- 
rium. We integrate over all radii to obtain the total de- 
cay rate idecay and, once again assuming steady state, 
take the star formation rate Mgp to estimate E- m from 
Eq. ([7]). This allows us to calculate the minimum efficiency 
e = |idccay |/iin needed for sustaining the observed disk 
turbulence by gas accretion. The results are presented in 
Table [3] and graphically illustrated in Figure [SJ 

We notice that all galaxies similar to the Milky Way, 

1. e. those with Hubble types ranging from Sb to Sc with 
v T ot 200 kms -1 , require efficiencies of only 10% or less. 
This holds despite variations in star formation rate or to- 
tal gas mass of almost a factor of ten. As in Section 
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Fig. 3. Radial distribution of 3D velocity dispersion a (in units of kms -1 ), HI column density Ehi {Mq yr" 1 ), cor- 
responding turbulent scale Ld (kpc) taken as twice the disk scale height, computed from hydrostatic balance using 
Eq. (|12[) (solid lines) and from the total potential using Eq. (fTB")) (dashed line), number of pixels as function of radius 
in annuli of width AR = 0.5 kpc in the 11 THINGS maps. The last column provides an estimate of the statistical 
significance of the data at each radius. The radius is scaled to the optical radius i?25- The thin lines indicate the 95% 
variation. 
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Fig. 4. Local kinetic energy per unit surface area -Ekin (crgpc~ 2 ), and local turbulent energy decay rate per unit 
surface area .Edecay (ergs -1 pc -2 ) using Ld derived from vertical hydrostatic balance, Eq. (|12p. together with the 
cumulative total kinetic energy and cumulative turbulent dissipation rate as function of normalized radius for 11 
THINGS galaxies. The thin lines in columns 1 and 2 indicate the 95% variation at each radius using the scale height 
derived from Eq. (|12p. 



10 Klessen & Hennebelle: Accretion-Driven Turbulence 

Table 2. Observed properties of the analyzed THINGS galaxies. 





Hubble 


D 


^rot 


i?25 


Mm 




(a) 


(Shi) 




(Ld)pot 




type 


Mpc 


kms - 


kpc 


10 9 M 


Msyr" 1 


kms 


M pc- 2 


kpc 


kpc 


Name 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V 


(8) 


(9) 


(10) 


NGC 5194 


Sbc 


8.0 


220 


9.0 


2.5 


6.05 


31.69 


7.18 


4.72 


1.32 


NGC 628 


Sc 


7.3 


220 


10.4 


3.8 


1.21 


14.28 


4.55 


1.00 


0.80 


NGC 3184 


Sc 


11.1 


210 


11.9 


3.1 


1.43 


18.34 


4.48 


1.46 


0.99 


NGC 3351 


SI) 


10.1 


200 


10.6 


1.2 


0.71 


20.89 


2.05 


5.42 


1.01 


NGC 6946 


Sc 


5.9 


200 


9.8 


4.2 


4.76 


18.72 


5.57 


1.28 


1.13 


NGC 5055 


Sbc 


10.1 


200 


17.2 


9.1 


2.42 


23.58 


3.15 


4.15 


2.81 


NGC 4736 


Sab 


4.7 


160 


5.3 


0.4 


0.43 


24.98 


2.99 


5.11 


0.71 


NGC 7793 


Scd 


3.9 


130 


6.0 


0.9 


0.51 


19.64 


6.06 


1.40 


0.91 


IC 2574 


Sm 


4.0 


80 


7.5 


1.5 


0.12 


17.30 


5.82 


1.13 


2.42 


NGC 4214 


Irr 


2.9 


60 


2.9 


0.4 


0.05 


16.91 


4.78 


1.18 


1.17 


HO II 


Irr 


3.4 


40 


3.3 


0.6 


0.07 


16.89 


6.37 


0.96 


1.99 



(1 ) H ubble type as listed in the LEDA data base (URL: |http://led a.uni v-lyonl.fr/). (2) Di stance according to IWalter et alj 
(|2008r ). (3) Peak of the rotation curve as obtained from the Appendi x in lde Blok et al.l (|2008h . Note that our values agree well 
with the fit formula provided in Appendix B.l in iLerov et al.l (|2008l ) for galaxies where the flat part of the rotation curve is 
observed. The only exception is IC 2574, which shows a continuously rising roation curve within the observed radius range. (4) 
R25 is the B-band isophotal radius at 25 mag arcsec -2 , which is a standard proxy for t he optical ra d ius of the galaxy. (5) 
Total HI mass of the galaxy, obtained by integrating ov er all pixels, see also Table 4 of ILerov etaH (l2008h . (6) Total star 
formation rate of the galaxy, as provided by Table 1 of IWalter et all (|2008l ) or Table 4 of ILerov et all ||2008h . (7) Surface- 
density weighted 3-dimensional mean velocity dispersion in the galaxy. (8) Mean HI surface density. (9) Maximum length scale 
of the turbulent velocity field. (Ld)HD is calculated from the disk thickness in each pixel as Ld = 2H with H = a 2 D / (27tG£ht) 
and then averaged over the entire galaxy using surface-density weighting. Note, that Ld varies strongly with radius (see Figure 
[3|, so that the physical interpretation of average turbulent length scale is not straight forward. (10) Estimate of the mean 
turbulent dissipation scale (-Ld)pot based on the disk thickness derived from the potential method, Eq. H13|) . Note that both 
approximations give roughly the same numbers with (L d )HD > (L d }pot for galaxies with w rot > 200 kms -1 and (L d )HD < (Ld)pot 
for the dwarfs with v rot < 100 kms -1 . 



Table 3. Derived energy input and decay rates and corresponding minimum efficiencies. 





Ein 


Ldccay.HD 


£hd 


-C/decay ,pot 


£pot 


E decay ,R.25 




Skin 


f -^kin 


/■^'dccay,HD 


e >«25 




10 39 ergs -1 


10 39 ergs -1 




10 39 ergs" 1 




10 39 ergs" 1 




10 54 erg 




Name 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V 


(8) 


(9) 


(10) 


(11) 


NGC 5194 


92.39 


3.41 


0.037 


7.08 


0.077 


0.45 


0.005 


28.20 


0.26 


0.17 


0.006 


NGC 628 


18.48 


1.97 


0.106 


2.06 


0.111 


0.06 


0.003 


8.39 


0.38 


0.28 


0.027 


NGC 3184 


19.90 


0.82 


0.041 


0.96 


0.048 


0.04 


0.002 


11.03 


0.21 


0.13 


0.005 


NGC 3351 


8.96 


0.16 


0.018 


1.08 


0.121 


0.06 


0.006 


7.12 


0.15 


0.17 


0.003 


NGC 6946 


60.07 


4.44 


0.074 


5.63 


0.094 


0.28 


0.005 


17.01 


0.32 


0.29 


0.020 


NGC 5055 


30.54 


1.56 


0.051 


3.59 


0.117 


0.15 


0.005 


58.98 


0.31 


0.15 


0.008 


NGC 4736 


3.47 


0.52 


0.151 


6.42 


1.848 


0.29 


0.083 


3.52 


0.17 


0.12 


0.018 


NGC 7793 


2.72 


1.87 


0.686 


1.95 


0.716 


0.13 


0.050 


3.57 


0.17 


0.07 


0.047 


IC 2574 


0.24 


2.18 


9.009 


0.60 


2.459 


0.10 


0.417 


4.52 


0.28 


0.18 


1.600 


NGC 4214 


0.06 


1.28 


22.457 


1.17 


20.549 


0.21 


3.653 


1.33 


0.49 


0.42 


7.831 


HO II 


0.04 


1.77 


50.070 


0.67 


19.065 


0.17 


4.900 


1.81 


0.59 


0.50 


18.554 



(1) Total kinetic energy provided by infalling gas, calculated from Eq. ^ using M = Msf and Vi B = v I0 t- (2) Total dissipation 
rate of turbulent kinetic energy L7decay,HD obtained from integrating Eq. ((2| over the entire galaxy using Ld derived from 
vertical hydrostatic equilibrium. (3) Minimum efficiency ehd = |-Edecay,HD|/-Ein required for the conversion of infall motion into 
turbulent energy in the disk. (4) Integrated turbulent decay rate -Edecay.pot based on the potential method. (5) Corresponding 
minimum efficiency e po t- (6) Integrated turbulent decay rate -Ed eC ay,R 2 5 based on the assumption that the outer scale of the 
turbulent velocity field is equal to the size of the disk Ld ~ 2 7?25- We list this value to provide an estimate of the uncertainty 
introduced by not knowing Ld very well. In Section l3~5"l we argue that -Edecay,R 25 is a strict lower limit and that the true decay rate 
probably lies more closely to -Edecay.HD- (7) The corresponding required minimum efficiency £r 25 for accretion driven turbulence 
to work. (8) Total kinetic energy Lki n integrated over the whole galaxy. (9) Fraction of total kinetic energy outside the optical 
radius i?25- (10) Fraction of turbulent energy decay rate outside of R25 using the disk vertical scale height from hydrostatic 
balance. The total rate is given in column (2). (11) Corresponding required accretion efficiency to drive the turbulence in the 
outer disk for R > i?25- Compare to the total galactic value given in column (3). The outer disk value is typically by a factor 
of 5 lower. 
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Fig. 5. Minimum efficiency required to sustain the ob- 
served level of turbulence by accretion from the galactic 
halo for the sample of THINGS galaxies. Our fiducial val- 
ues, £hd (filled circles) are based on using the disk thick- 
ness from vertical hydrostatic balance as outer scale of 
the turbulent velocity field To give an illustration of 
the uncertainties involved in estimating Ld, we also show 
e P ot derived from using the global potential to calculate 
the disk thickness (horizontal symbol), and £r 25 where we 
take the diameter of the optical disk as proxy for (ar- 
row down). The latter is a strict lower limit, with the true 
minimum efficiencies probably being closer to our fiducial 
estimate. Note, that for the dwarf galaxies in our sample 
(vmt < 100 km s" 1 ) ejjD ^ £pot while the opposite is true 
for the Milky Way type galaxies (v TO t <; 200 km s -1 ). 



we argue that e < 0.1 is a reasonable number. Given that 
the true turbulence dissipation scale could exceed the disk 
thickness and potentially be as large as the disk diameter 
(see Section 13.51 below) , the derived decay rates are up- 
per limits and the minimum efficiency of accretion driven 
turbulence might be considerably smaller. We point out 
that the model clearly fails for the dwarf galaxies in our 
sample with w rot < 100 km s" 1 . The kinetic energy added 
by accretion is not sufficient to compensate for the energy 
loss by turbulent decay. Either these galaxies accrete more 
mass than inferred from their low star-formation rates, 
or there are other processes that dominate the disk tur- 
bulence on all scales. Both explanations appear equally 
likely and ask for more detailed studies with specific focus 
on dwarf galaxies. 

In order to better understand the physical origin of this 
variations, we plot the minimum efficiency values obtained 
above as function of different galaxy parameters in Figure 
[H We see that e correlates very well with the rotational 
velocity v rot of the galaxy, log 10 e ~ 2.3 — u ro t / (60 km s - 1 ) . 
This is understandable because we use v ro t as proxy for 
the infall velocity and hence for the kinetic energy of the 
infalling material. The rotational velocity is a good mea- 
sure of the total mass including dark matter and stellar 
component. If we exclude the dwarf galaxies in the sam- 
ple, we find no correlation between e and total HI mass, 



star formation rate, and velocity dispersion. This is some- 
what surprising. We had expected the efficiency to scale 
inversely with HI mass, because the total amount of tur- 
bulent energy that needs to be replenished scales linearly 
with Mhi assuming that a is roughly constant. The con- 
trary is the case, if we focus on the dwarf galaxies we see 
that they have high e values despite small Mm. Because 
the decay rate, Eq. (|2|), scales with the third power of the 
velocity dispersion, we had also expected that e strongly 
depends on a. Again, this is not seen. From our small 
sample of 11 galaxies, we conclude that the total mass of 
the galaxy is the main parameter determining the poten- 
tial importance of accretion driven turbulence. It would 
be interesting to perform a similar analysis with a larger 
sample of galaxies to have better statistics. 

3.5. Main Uncertainties 

The processes discussed here are subject to large uncer- 
tainties. First of all, virtually all quantities that enter the 
theory vary with radius. Defining a galaxy- wide average 
value is not a trivial task. In cases where we have extended 
HI maps, this is not a problem, because we can integrate 
over all radii and obtain well defined global values for E^, 
-E'decay, an d E m . However, for galaxies with only a few pix- 
els across or for observations with insufficient sensitivity to 
detect the extended disk the errors can be considerable. 
Second, we do not know how the gas enters the galaxy. 
Docs it fall onto the outer disk in discrete cold streams 
(as indicated by cosmological simula tions at high redshift 
Dekel et al.ll2009HAgertz et al.ll2009t) ? Or does it condense 
out of the tenuous halo gas and enter t he disk more gently 
and more distributed as proposed by Peek! ( 20091 )? Both 
processes could lead to very different efficiencies. 

The third and probably most severe uncertainty con- 
cerns the turbulent length scale. The observational data 
are not very conclusive, with estimates of the turbulent 
scale ranging from only 4pc ( Minter fc Spangler 1996t) 
from s cintillation measurem ents in the Milky Way up to 
6 kpc (|Dib fc Burkertll2005f ) from computing the autocor- 
relation length of HI in Holmbcrg II with large uncer- 
tainties in both values. Very careful statistical analyses 
of the power spectrum in se veral nearb y molecular c louds 
bv lOssenkopf fe Low! (|2002l ) as well as iBruntl (|2003l) and 
Brunt et all (|2009h . however, indicate that the bulk of the 
turbulent energy is always carried by the largest scales 
observed. This is consistent with turbulence being driven 
from the outside. 

Throughout most of this paper we assume that the 
outer scale of the turbulent cascade is comparable to 
the disk thickness, i.e. we take it as being twice the 
vertical scale height. We argue that this is the relevant 
upper length scale in the system. Only then can we 
speak of 3-dimensional turbulence, whe re we have good 
estimates of the decay properties (e.. 



Stone et al 
ElmegreerT 



19981 IPadoan fc Nordlun 



.Low et al 
19991: 



Low 



1998 



1999; 



2000). One may propose, however, that the 
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Fig. 6. Correlation between required minimum efficiency e and rotational velocity v TO t (top left), total HI mass Mm 
(top right), star formation rate (bottom left), and average 3-dimcnsional velocity dispersion a (top right). The used 
symbols are the same as in Figure [51 Note that in the top left plot, we use open circles to denote the shift introduced 
when using the fit parameter i>fl at 



Lerov et al. I (120081) instead of the observed maximum of the rotation curve (Table 



tab:THINGS). This is only significant for the galaxy IC 2574 which has a continuously rising roation curve within the 
observed radius range. The efficiency e correlates well only with v mt . 



turbulent cascade extends all the way across the galac- 
tic disk. In this case, ~ 2i? 2 5- Turbulence is mostly 
2-dimensional and it is not clear how to estimate its de- 
cay properties in the differentially rotation disk. It has 
been speculated, however, that the decay properties of 
2- and 3-dimensional turbulent flows could be equivalent 
as long as they are in the strongly supersonic regime, 
because dissipation occurs mostly in s heet-like shocks 
(jAvila-Reese fc Vazouez-Semadenj l200lh . If Eq. © is 



still approximately correct, the decay rate is considerably 
lower because R25 ^> H, and our model can tolerate even 
lower minimum efficiencies (Table [3]). On the other hand, 
the velocity difference across the disk is of order of the ro- 
tational velocity v ro t , which for the Milky Way type galax- 
ies greatly exceeds a. Because H/R ~ cr/v Ioi , Eq. (TT3| . it 
is likely that both effects approximately cancel and that 
the decay rates based on H and on R25 are comparable. 
This is in fact what we would expect from a self-similar 
turbulent cascade, where the behavior is determined by 
the physical properties at the dissipation scale. We adopt 
Ld based on the disk scale height obtained from vertical 



hydrostatic balance, Eq. (TT2"j) . as our fiducial value. The 
range of e values for each galaxy in our sample associated 
with the uncertainties in is illustrated in Figure [5] 

We also want to call attention again to the fact that we 
have neglected the molecular gas content in our analysis of 
the THINGS galaxies. Indeed, some of these galaxies con- 



tain an appreciable amount of molecular gas ([Lerov et al 



20091) . However, this gas is mostly contained within R25 
and in addition carries only little turbulent kinetic energy 
compared to the atomic component (see also Section 13.31 
for the Milky Way). The error involved in focusing on 
atomic gas only is therefore small. 

3.6. Outer Disks of Spiral Galaxies 

It is well know now that the g aseous disk extends 
well beyond the optical radius ([Thilker et al.l 12007 : 
Zaritskv k Christleinl 120071) . The sample of THINGS 



galaxies allows us to study the turbulent energy content 
of extended disks in more detail. We begin by asking 
what parts of the disk carry most of the turbulent ki- 
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netic energy. We calculate the cumulative kinetic energy as 
function of radius as well as the cumulative energy decay 
rate and plot both quantities columns 3 and 4 of Figure 
I3J We see that a significant fraction (between 15% and 
59%) of the total kinetic energy is carried by the outer 
disk. The numbers together with the required accretion 
efficiency are provided in Table [3] This finding has impor- 
tant consequences for our understanding of the origin of 
this turbulence. While within R25 energy and momentum 
input from stellar sources (supernovae, stellar winds and 
outflows, expanding HII regions) can contribute signifi- 
cantly to driving ISM turbulence, this approach fails for 
the outer parts. Here accretion from the extended gaseous 
halo, ma ybe in concert with t he magneto-rotational in- 
stability ( Tamburro et al. 20091 ). is the only astrophysical 
driving source available. The required efficiencies for Milky 
Way type galaxies are about 1%. This is a very low value 
and we conclude that accretion could easily drive the tur- 
bulence in the outer disk of present-day spiral galaxies. 



l/2vf n w GM/L d and Eq. © simplifies to 



1^ 

e 2G 



(15) 



With typical clump sizes Ld ~ 1 kpc, masses of about 
~ 10 8 M Q , and 3-dimensional velocity dispersions a ps 
30kms , the required accretion rate to drive the clump's 
internal turbulence is 



M in w - 1 M yr" 1 



(16) 



The observed star formation rates of clumpy galaxies lie 
in the range 10 — 50M Q yr~ 1 . If we assume that these 
numbers can be used as proxy for the mass accretion rate 
then we see that again efficiencies of order of 10% or less 
are sufficient to drive the clumps internal turbulence. This 
is consistent with the e-values estimated from numerical 
experiments as discussed in Section l2~2l and adds further 
support for the hypothesis of accretion driven turbulenc e 
in clumpy galaxies (see also Elmegreen fc Burkert 20091 ) . 



3. 7. Clumpy Galaxies at High Redshifts 

The mechanism that we propose here to drive internal 
turbulence is very generic and likely to operate on many 
different spatial and temporal scales. High rcdshift galax- 
ies, for example such as detected in the Hubble Ultra Deep 
field, are observed to be very irregular, with blue clumpy 
structure, asymmetry, and a lack of cen t ral concentration 
dElmegreen et al.l 120051 : IConselice 2003 : lElmeereen et al 



2009b). They are typically characterized by a consider 
ably higher degree of internal turbulence, e.g. as reflected 
by the large observed line width of Ha e mission, than 



prese nt-day galaxies in the same mass range ([Genzel et al 



20081 ). There seems to be an evolutionary trend with de- 
creasing rcdshift from clumpy galaxies with no evidence 
of intcrclump emission to those with faint red disks. This 
trend continues to prese nt-day spiral galaxies o f occu- 
lent or grand design types (jElmegreen et al.ll2009bl ) . Some 
clumpy galaxies at high redshift resemble massive versions 
of local dwarf irregular galaxies. They exhibit very high 
gas fractions and appear to be in a very young evolu- 
tionary state. Soon after their discovery it has been rec- 
ognized that their strong turbulence is difficult to main- 
tain with internal sources as the stellar feedback processes 
that act in present-day galaxies are relatively ineffec- 
tive when the vel ocity dispersion of the w hole interstellar 
medium is large dElmegreen et al. Il2009ah. Instead it has 
been p roposed (e.g. Genzel et al.l 20081 : Elmegreen et al.l 
2009bJ) that this high degree of turbulence could be 



driven by cold accretion streams as found in numeri- 
cal simulations with detailed treatm ent of the thermo- 
dynamic behavior of the infalling gas ([Birnboim fc Dckcl 
20031: ISemelin fc Combed 120051: iDekel fc Birnboim! 120061: 
Agertz et al.ll2009t ICeverino et al.|]2009h . 



3.8. Further Discussion 

The analysis above relies on the assumption that galax- 
ies evolve in steady state so that the star-formation rate 
is matched by the infall of fresh gas from an extended 
halo. For most present-day field galaxies this appears to 
be a reasonable assumption. However, it breaks down for 
highly perturbed systems, e.g. when galaxies experience a 
major merger or are tidally disturbed in the central regions 
of dense galaxy clusters. Strong perturbations lead to en- 
hanced star formation without being necessarily accompa- 
nied by additional gas infall. It is the original disk gas that 
is converted into new stars at increased rate. We therefore 
expect our model to work best for disk galaxies that are 
marginally unstable and where stellar birth proceeds in 
a self-regulated fashion similar to the Milky Way. Indeed 
the L± galaxies in our sample with w rot ~ 200 km s -1 all 
fit our model very well, while the dwarf galaxies with 
Vrot ^ 100 km s -1 do not. These galaxies are character- 
ized by very irregular clumpy structur e. They resemble 



the c l umpy galaxies se en at high redshift (jElmegreen et al 



20051: IConselicell2003h . Taken at face value, our simple 
equilibrium model therefore should not apply to those 
galaxies either, as they most likely are still in the phase 
of rapid mass growt h through merging and accretion of 
massive cold clouds ( Genzel et al . 2008: lElmegreen et al.l 
2009^). However, we can apply our model to structures 



If we assume the accretion onto the clump at any in- 
stance in time is driven by the clump's self-gravity, then 



within these galaxies and speculate that the turbulence 
observed within the dense clumps is driven by accretion 
fSection !3.7p . This is similar to the mechanism we propose 
in Section U to drive turbulence in present-day molecular 
clouds in the Milky Way. 

We also need to point out, that our analysis by no 
means implies that other sources of turbulence are not 
important. Clearly in the inner parts of the disks of L+- 
typc galaxies, stellar feedback plays a key role and can 
provide sufficient energy to drive the observed turbulence 
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(jLow fe Klessenll2004l ). ] n addition, it is difficult to see how 
the accretion energy, that we expect due to angular mo- 
mentum conservation to be mostly added to the outer re- 
gions of the disk, is able to reach the inner disk within rea- 
sonable timcscalcs and without being dissipated. If asso- 
ciated with net mass transport, one would expect a mean 
inward flow of order of 5kms - 1 , which is not observed 
(Wong et al. 2004, however, see |Peekll2009l for an alter- 
native accretion scenario). It is important to note further, 
that the accretion of halo gas is likely to be non-spherical 
and may excite wave-like perturbations that then could 
tap to the system's rotational energy. This provides an 
additional reservoir for driving turbulence in the disk so 
that even small accretion rates could lead to effective disk 
heating. 

Furthermore, we note that our analysis neglects the 
influence of galactic fountain flows. Expanding supernova 
bubbles or large HII regions could transport hot and metal 
enriched material into the halo, where it cools and eventu- 



ally f alls back onto the galactic disk (jCorbelli fe Salpeter 
19881) . This process is a key element of the matter 
cycle and enrichme nt history on global galactic scales 



( Spitoni et al. 2009). In addition, the kinetic energy as- 



sociated with the returning material could contribute to 
driving ISM turbulence just like infalling fresh material 
does. This supplementary source again renders the above 
estimate of the minimum efficiency required for accretion- 
driven turbulence to work an upper limit. In particular, 
galactic fountains should not affect the outer disk much, 
because they deliver the ejected material close to the ra- 
dius where it originates from. They are i mportant only 



of turbulence in molecular clouds with no, or extremely 
low star formation like the Maddalena cloud or the Pipe 
nebula, is significant and broadly comparable to the level 
of turbulence observed towards star forming clouds. Both 
facts seem difficult to reconcile with turbulence being 
driven from internal stellar sources. 



We argue that it is the very process of cloud formation 
that drives its internal motions by setting up a turbulent 
cascade that transports kinetic energy from l arge to small 
scale s in a universal and self-similar fashion ([Kolmogorov 
19411 ). Our hypothesis is that molecular clouds form at 
the stagnation points of large-sc a le convergent flows (e.g . 
Ballesteros-Paredes et al.l 19991 Hartmann et al. 2001 



for the inne r, star-forming parts of the disk (jMelioli et al 
2008ll2009l) . 



Klessen et al.ll2005tlHeitsch et al.ll2006bHfIennebelfe et al 
20081 : iBaneriee et aDl2009l ). maybe triggered by spiral den 
sity waves or other global perturbations of the gravita- 
tional potential. As the density goes up the gas can cool 
efficiently, turn from being mostly atomic to molecular, 
and shield itself from the external radiation field. As long 
as the convergent flow continues to deliver fresh material 
the cloud grows in mass and is confined by the combined 
thermal and ram pressure of the infalling gas. Because 
the molecular gas is cold its internal turbulent motions 
are strongly supersonic. Consequently, the cloud develops 
a highly complex morphological and kinematic structure 
with large density contrasts. Some of the high-density re- 
gions become gravitational unstable and go into collapse 
to form stars. This modifies the subsequent evolution, as 
stellar feedback processes now contribute to the energy 
budget of the cloud. 



4. Molecular Cloud Turbulence 

4.1. Theoretical Considerations 

Since the seminal work of LarsonI ( 198ll) . it is well es- 
tablished that molecular clouds are highly turbulent. The 
3-dimensional velocity dispersion in these objects varies 
with their size, L, and typically follows the relation, 



a 0.8 kms 



( 



L 



0..i 



(17) 



The physical origin of this turbulence is not fully 
understood yet. In particular, the question as to 

whether it is i njected from the outside, e.g. by 

colliding flows ( Hunter et al. 19861 Vishniac 1994 ). 
or driven by internal sources such as protostellar 
outflows (ILi fe Nakamural l2006t IBaneriee et al.l 12007 



Nakamura fe Lill2008HWang et al.l |2010|) or exp anding HII 
egio ns ( MatznerT 2002h or supernovae ( Low fe Klessen 



20041) , is still subject to considerable debate. We favor the 
first assumption as observations indicate that molecular 
cloud turbulence is always domin ated by the largest-scale 
modes accessible to the telescop e (jOssenkopf fe Lowll2002l: 
Bruntll2003t iBrunt et all 120091 ) . In addition, the amount 



A long series of numerical simulations focusing on 
molecular cloud dynamics and attempting to build 
up these clouds from diffuse gas at the stagnation 
points of convergent larger-scale flows have indeed shown 
that accretion can sustain a substanti al degree of tur 



bulence in the newly formed cloud (Waldcr & Folini 



1998 1 120001: iKovama fc Inutsukal l2002t 



1998 




2006 


: 1 


2006 


1 - 



Folini fc Waldei 



2006 ; Heitsch et al.ll2006at IVazquez-Semadeni et al.ll2003 



Figure [7] shows the internal velocity dispersion of molec- 
ular cloud clumps extracted from the high resolution nu- 
merical simulation described in Section 12.21 as a function 
of their size. These clumps are defined by a simple clip- 
ping algorithms, selecting cells with densities larger than 
2500 cm -3 . The internal velocity dispersion is then com- 
puted by computing the rms velocity with respect to the 
cloud bulk velocity. As can be seen, the velocity dispersion 
is compatible with Larson's relation, Eq. (fTT)) . To further 
illustrate this, Fig. [5] shows the column density of one of 
the clumps formed in the simulation. These simulations 
suggest that continuous accretion of diffuse material in a 
molecular cloud is sufficient to maintain a high level of 
turbulence inside the cloud. 
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Fig. 7. Internal velocity dispersion of clumps produced in 
colliding flows simulations. 




Fig. 8. Column density for a clump produced in the sim- 
ulation. 



4.2. Application to Molecular Clouds in the LMC 

The best evidence for gas accretion inside molecular clouds 
may be given by observa tions in th e Larg e Magellanic 



Cloud fLM C) reported bvlBlitz et all (|2007l) . iFukui et al 
(|2009h . and lKawamura et al.l ( 2009 ). These authors distin- 



guish 3 types of giant molecular clouds that they interpret 
as an evolutionary sequence. During the first phase, which 
should last about 6 Myr based on statistical counting, the 
clouds are not forming massive stars and therefore exhibit 
a low star formation rate. During the second phase, which 
may last about 13 Myr, massive stars form but not clus- 
ters. The last phase is characterized by the presence of 
both massive stars and clusters. While the mean mass of 
the clouds observed in the first phase is of the order of 
1.1 x 10 5 Mq the mean mass of the clouds in the second 
phase is about 1.7 x 10 5 M Q . This implies that the giant 
molecular clouds in the LMC are on average accreting at 
a rate of about 



M « 1O" 2 M yr" 



(18) 



Note that Fukui et al. ( 2009h quote a slightly larger value 
of about 5 x 1CP 2 M Q yr -1 and note also that similar 
values are found in numerical simulations of molecular 
cloud formation in convergent flows (see, e.g. Figure 7 in 
Vazquez-Semadeni et al. I l2009h . In the following discus- 



sion, we adopt our estimate as fiducial value, but keep in 
mind the alternative higher estimate. 

The amount of energy which is delivered by this pro- 
cess as well as the energy dissipated per unit time in the 
cloud are given by Eq.'s. © and To obtain an esti- 
mate for the relevant parameters, we take the outer scale 
of the turbulence, L^, to be equal to the size of the molec- 
ular cloud, and adopt the observed relation between clou d 
mass and size (see, e.g., Figure 1 in lFakarone et al.ll2004l ). 



M « 10 M 



Vlpc 



2.3 



(19) 



Note, that this empirical behavior implies that the cloud 
has a fractal dimension of 2.3 or, more or less equivalently, 
has large internal density contrasts. Recall that the scal- 
ing relation for homogeneous clouds is simply M oc L 3 . We 
then calculate the infall velocity V[ n using mass conserva- 
tion. The accretion rate M must be equal to the flux of 
mass through the cloud surface. Assuming spherical sym- 
metry this surface is 4irR 2 leading to the relation 



M = 4irR UjnPisM, 



(20) 



where Pism is the typical ISM density outside the cloud. 
Note that the assumption of spherical symmetry gives 
an upper limit for v- ln since it minimizes the surface. 
Consequently, our estimate of E m is also an upper limit. 
With Eq. we obtain 



15 km s 



M 



10" 2 M yr- 



. cm 



-3 



La 



80 pc 



(21) 



where is scaled to the diameter of a molecular cloud 
with mass ~ 2 x 10 5 Mq and where Hism is the number 
density of the atomic gas using a mean molecular weight 
of 1.3. Thus, it appears that flows with velocities of about 
15kms _1 are sufficient to explain the observed accretion 
rate of about 10~ 2 M Q s _1 . 

Combining Eq.'s @ and §5§ we obtain 



E, 



e = 



decay 



7r 2 pj SM ML"^ M 



3^3 



which translates into 

2 







04 (™Y 

VI cm A J 



Al 



10 5 M m 



M 



10- 2 Meyi- 1 



(22) 



-3 

(23) 



with the help of Eq.'s (IT) and (fl9"]l. 

The kinetic energy associated with formation and sub- 
sequent growth of molecular clouds is sufficient to drive 
their internal turbulence provided that the efficiency of 
this conversion is not smaller than a few per cent. These 
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numbers are in very good agreement with the numerical 
results discussed in Section 12.21 with mean density con- 
trasts between cloud and intercloud medium of several 
10s to 100. Taking into accou nt the higher accretion rate 
quoted by iFukui et al. allows for even small effi- 

ciencies or equivalently for larger density contrasts. 

5. Turbulence in Accretion Disks 

Last we investigate whether accretion onto T Tauri disks 
may represent a significant contribution to the turbu- 
lence in these objects. We focus our discussion on the 
late stages of protostellar disk evolution, the class 2 and 
3 phase, where the original protostellar core is almost 
completely accreted onto the central star, which is then 
surrounded by a remnant disk carrying only a few per 
cent of the total mass. Because the system is no longer 
deeply embedded, it is accessible to high-precision obser- 
vations and structu re and kinematics are well constrained 
(jAndre et al.ll2000l) . Although it is currently thought that 
turbulence in protostellar accretion disks is driven through 
the no n- linear evolution of th e magneto-rotational insta- 
bility (|Balbus fe Hawlevlll998l ). we believe that it is never- 
theless worth to estimate the level of turbulence that the 
forcing due to accretion may sustain. 

5.1. Are T Tauri Disks Accreting ? 

The accretion of gas onto the disk is not easy to measure 
during the T Tauri phase and no observational data are 
available in the literature. On the other hand, the accre- 
tion from the disk onto the star has been measured in a va- 
riety of objects. Typical accretion rates are of the order of 
3 x 10~ 8 M ff yr -1 for a 1 solar mass star (e.g. iNatta et al 



2004 IMuzerolle et~alll2005l: lLopez et al.ll2006t lGatti et al 



2006, 2008) while the mass of the disk is typically of the 
order of 10 -2 M Q . Dividing the latter by the former, we 
find that the disk could not last more than 3 x 10 5 yr 
whic h appears to be sh orter than the typical T Tauri ages 
(e.g. lEvans et al.ll2009l) . Thus, it seems likely that T Tauri 
disks are still accreting gas at a rate comparable to the 
one at which gas fr om the disk is accr e ted onto the cen 
tral star 



( sc ^ e.g. 



Padoan et al. 2005; Dullemond et al 



20061: iThroop fe Ballvll2008h . Alternative ly, the disk could 



be m ore massive than usually assumed (IHartmann et al 
20061) . 



Another interesting observation is that the accretion 
rate onto the star, M, is typically related to the stellar 
mass as 



M oc M" , 



(24) 



where a rss 1.8, howev er, with a large scatter of about one 



order of magnitude (IMuzerolle et al.l 2003 : Natta et al 



20061) . IPadoan et al 




2005D and IThroop fe Ballvl (|20~ 

propose that the stars as they move through the 
cloud, accumulate ga s at the Bondi-Hoyle accretion rate 
( Bondi fc Hovlejl944 ). Estimating the velocity of the star 



to be about 1-2 km s" 1 and taking a mean cloud den- 
sity of 10 3 " 4 cm" 3 , they infer typical accretion onto a one 
solar mass star of about 10 -8 Mqvt . In particular, to 
get accretion rates compatible with the largest observed 
values, they need to invoke large densities of 10 4-5 cm" 3 . 
However, the filling factor of gas at such densities is very 
low, of order of a few per cent. Densities of 10 4-5 cm" 3 
typically correspond to cores with sizes of ~ 0.1 pc. If we 
take typical stellar velocities of 1 — 2kms -1 , it requires 
only 10 5 yr to cross a core. Most of the time T Tauri stars 
travel through low density gas, eventu ally building up a 



widely dispersed population (see, e.g. iNeuhaeuser et al 



1995t Iwichmann et all fl 9961) 



Here we propose that accretion during the late phases 
of protostellar evolution proceeds in a different way. We 
base our discussion on the assumption that the velocity 
of the star is inherited from the bulk velocity of the core 
in which it forms. In p-Oph, for example, the typical 
core-to-cor e velocity dispersio n is found to be less than 
0.4 km s" 1 (I Andre et alj|2007h . The typical stellar veloc- 



ity dispersion measured in nearby T Tauri associations 
or open star clu sters is of sim ilar order, ~ 0.3 km s" 1 
for the Hyades (jMadsenl 120031) . ~ 0.5 - 0.6 km s" 1 for 
Com a Berenices, Pleiades, and Praesep e dMadsen et al 
20021) . and below ~ lkms' 1 for a Per (jMakarovl 120061 ) 



Lupus jM akarov 120071 ) and the sub-groups in Taurus 
(|jones fc Herbidll979t feertout fc Genoval 120061) . It gets 
above lkms" 1 o nly for the more m assive clusters and 
OB associations (jMadsen et al.1 120021 ). Since star form- 
ing cores are part of a turbulent molecular cloud, we 
resort once again to Larson's relation and assume the 
mean velocity of any fluid element with respect to the 
center of the core is increasing with distance from the 
core as ~ 0.46 x (i/lpc) 5 . This expression is identi- 
cal to Eq. (JT7J) except for the factor 3" 1 / 2 which comes 
from the fact that we consider ctid instead of a. It takes 
a star about 2.2 x 10 6 yr to reach a distance of 1 pc if it 
travels with 0.46kms _1 . This implies that during a long 
period of time, comparable to the age of the T Tauri 
star, its velocity with respect to the surrounding gas is 
not of the order of 1 — 2kms~ 1 but more comparable to 
the sound speed of the gas c s ~ 0.2 km s" 1 . As we show 
below, we can quantitatively reproduce the observational 
relation M ~ 2 - 3 x 10~ 8 Moyr" 1 (Af*/1 M Q ) 18 ob- 
tained bv lNatta et al.l ( 20061) if we assume that the accre- 
tion onto the star is essentially controlled by the accretion 
onto the disk from the turbulent cloud environment. In 
this picture, the disk only acts as buffer for this overall 
acc retion flow (possibly leadin g to occasional outbursts, 



see 



Hartmann fc Kenvonlll996 ). 



Our esti mate is very c lose to the spherical accretion 
considered bv lBondil (1952). The difference however is that 
turbulenc e velocity increases with the distance (see also 
Rov 20071 for accretion in fractal media) . Consider a star 
of mass M* inside a turbulent cloud. Assuming that the 
star is at rest, it will be able to accrete gas inside a sphere 
of radius i? a cc, such that fluid particles inside this radius 
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are gravitationally bound to the star, 
^ GM+ 



(25) 

where c s w 0.2 kms is the sound speed and g\d = 
S" 1 / 2 ^ is the 1-dimensional velocity dispersion, which we 
obtain from Larson's relation, Eq. (|17p , assuming isotropy. 
From this relation we can easily obtain i? acc as a function 
of mass 




2vl 



(26) 



where vq m 0.8 kms -1 and Lq = lpc. The total mass 
included into the sphere of radius R aC c is simply M acc = 
47r/3i? acc p where p is the mean density inside molecular 
clouds, which we take to be of the order of /xxlOO cm~ 3 
with fj, = 2.35 x mfj being the mean molecular weight and 
mfj being the mass of the hydrogen atom. This mass falls 
into the star/disk system on a timescale of the order of 
Tacc ~ (GM+/ i? acc ) -1 / 2 leading to a typical accretion rate 
of 



W*&£ = (An/3)RU" c Ml /2 G^p. 



(27) 



Note that as the matter inside i? acc is accreted, the pres- 
sure drops and it is replenished in about a crossing time 
by the surrounding gas. Since by definition of i? acc , this 
crossing time is comparable to the accretion time, there 
is always fresh gas available for accretion. Altogether, we 
obtain 

Af acc oc Ml /2 (-1 + y/1 + £M*) , (28) 

with £ = 4Gv 2 /(3L a cj). 

In the limit where the cloud is not turbulent (a 2 <C c 2 , 
i.e. £ <C 1) this leads to Ma M 2 , while if it is dominated 
by turbulence (a 2 > c 2 or £ > 1) we obtain M oc M^ 4 . 
The slope of the M acc - M+ relation therefore lies between 
1.25 and 2. For a 1 solar mass star, R acc w 0.1 pc and 
a ps 0.25 kms -1 » c s , leading to a slope of ~ 1.8. Figurc[5] 



shows the dependency of M as specified by Eq. ([27)) (solid 
line). As molecular clouds are highly inhomogeneous the 
local density can vary over orders of magnitude. To ac- 
count for this variation, we also display the values of the 
accretion rate as predicted by Eq. l[2"7]) for n = 10 and 1000 
cm~ 3 (dashed lines). To connect to real measuremen ts, we 



overplot the data presented by lLopez et al. ( 20061) . The 
agreement is remarkable. Despite its simplicity, our model 
provides a good order-of-magnitudc fit. In particular, also 
our prediction that for masses between 1 and 10 M Q the 
relation between accretion rate and stellar mass becomes 
more shallow seems to be confirmed by the observational 
data. 

Finally, we stress that accretion of the type we consider 
here, is difficult to avoid as long as the star remains within 
its parent molecular cloud and we conclude that continu- 
ous accretion onto the disk even during the T Tauri phase 
is an interesting phenomenon that deserves further atten- 
tion. 




Kk-'' 

-0.5 0.0 

log(M,) [M B ] 



Fig. 9. Prediction of the accretion rate onto the disk 
as a function of the mass of the star. The solid line 
corresponds to a mean density of n = 100 cm~ 3 while 
the two dashed lines are for n = 1000 cm~ 3 (up- 
per curve) and n = 10 cm~ 3 (lower curve). To guide 
your eye the dotted lines indicate the slope of the re- 
lations M oc M 2 and M oc M+. We compa r e wit h 



data from ICalyet et all (120041). iMohantv et al l (|2005h 



iMuzerolle et all (120051) and iNatta et all (|200a as dis- 
played in Figure 3 of lLopez et al.l (|2006l ) , where crosses in- 
dicate detections and arro ws upper limits. The dot-dashed 
line is the fit proposed by INatta et al.l (|2006j ) with slope 
a = 1.8. 

5.2. Expected Velocity Dispersion in Accretion Disks 

Considering Eq.'s. © and ((3|), we again take the typical 
turbulent length scale to be comparable to the disk thick- 
ness, Ld = 2H, where H is the vertical scale height. Since 
the gravitational potential is dominated by the central 
star, we estimate the infall velocity v- m to be 



/ 2GM* 



1/2 



and get 



-E^dccay _ -^diskO' 3 /2iJ 



Mvi 



(29) 



(30) 



3 3/2 



M 6 



M 
4/G. 



M, 



disk 



where and Mdisk are the masses of the central star and 
its disk, respectively, and where A4 = ctid /c s is the rms 
Mach number of the turbulence in the disk. 

For a temperature of 10 K, c s rj 0.2 kms" 1 and the 
collapse rate c 3 /G w 2 x 10" 6 M o yr" 1 . With the fiducial 
value M « 2 x 10 -8 M yr _1 for a solar mass T Tauri star, 
the expression in the first bracket in Eq. (|30p is roughly 
100. In the class 2 and 3 phase, the ratio Mdi s k/M* « 0.01. 
We can estimate the ratio of local scale height and radius 
using Appendix |21 which gives H/R « 0.1 for M = 1 M Q 
and a disk radius of 200 AU. Together, this leads to 



24 M~ 



(31) 



1 /3 

Because A4 « (e/24) ' , late accretion can only drive sub- 
sonic turbulence in T Tauri disks. If the true efficiency 
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of the process is of the order of 0.1, then the maximum 
rms Mach number that can be reached is M. sa 0.16. If 
e = 0.01 then Ai w 0.07. Only few measurements of the 
velocity dispersion i n disks have been carried out to date. 
Dutrev et al.l (|2007l) report values of ~ 0.1 km s" 1 which 
corresponds to a Mach number M. sa 0.5. Although this 
value is slightly larger (factor 3 to 7) than the ones we 
inferred, there are large theoretical and observational un- 
certainties. We have not taken into account that accretion, 
particularly if it is non-axisymmetric may help to tap en- 
ergy from the rotational energy available in the disk. We 
conclude that our simple estimates support the hypothe- 
sis that accretion can provide at least some level of turbu- 
lence in circumstcllar T Tauri disks. Note that even if it 
can sustain a significant amount of turbulence, the ques- 
tion as to whether this turbulence could transport angular 
momentum efficiently is entirely open. 

Because M varies with M* as illustrated in Figure [9j 
the ratio Ed ecay / E- ln depends on M* as well. According to 
Eq. ([3"0"]l. the rms Mach number of the accretion-driven 
turbulence for any given efficiency is smaller in disks 
around low-mass stars than around high-mass stars. For 
a low-mass star with M = 0.1 an efficiency of 10% 
would correspond to M. = 0.02 only, while for M = 10 M Q 
e = 0.1 would lead to Ai = 0.2, assuming the standard 
values H/R = 0.1 and Mdi s k/M* = 0.01 in both cases. 

6. Conclusion 

When cosmic structures form, they grow in mass via accre- 
tion from their surrounding environment. This transport 
of material is associated with kinetic energy and provides 
a ubiquitous source of driving internal turbulence. In this 
paper we propose that the turbulence that is ubiquitously 
observed in astrophysical objects on all scales at least to 
some degree is driven by this accretion process. To sup- 
port our idea, we combined analytical arguments and re- 
sults from numerical simulations of converging flows to es- 
timate the level of turbulence that is provided, and applied 
this theory to galaxies, molecular clouds, and protostcllar 
disks. 

We first studied the Milky Way as well as 11 galaxies 
from the THINGS survey, and found that in Milky Way 
type galaxies the level of turbulence ubiquitously observed 
in the atomic gas in the disk can be explained by accretion, 
provided the galaxies accrete gas at a rate comparable to 
the rate at which they form stars. Typically, the efficiency 
required to convert infall motion into turbulence is of order 
of a few percent. This process is particularly relevant in the 
extended outer disks beyond the star-forming radius where 
stellar sources cannot provide alternative means of driving 
turbulence. It is attractive to speculate that the popula- 
tion of high-velocity clouds, e.g. as observed around the 
Milky Way, is the visible signpost for high-density peaks 
in this accretion flow. The assumption of steady state evo- 
lution, however, fails for dwarf galaxies. In order to drive 
the observed level of turbulence the accretion rate needs 
to exceed the star formation rate by far and we expect 



other sources to dominate. We also applied our theory 
to the dense star-forming knots in chain galaxies at high 
redshift and, in agreement with previous studies, came to 
the conclusion that their turbulence could be driven by 
accretion as well. 

We then turned to molecular clouds. Using the re- 
cent estimate by Fukui et al. (2009) of the accretion rate 
within molecular clouds in the large Magenallic cloud, 
we found that accretion is sufficient to drive their in- 
ternal turbulence at the observed level. This is in good 
agreement with the finding that most of the turbulent 
energy i n molecular clouds carried by large - scale modes 
see, e.g. lOssenkopf fc Lowll2002l : brunt 20031 iBrunt et al" 



2009h and also with the fact that clouds which do not 
form massive star show the same amount of turbulence as 
those which do. This excludes internal sources. It is the 
very process of cloud formation that drives turbulent mo- 
tions on small scales by establishing the turbulent cascade. 
Numerical simulations of colliding flows reveal that the 
turbulence within dense clumps generated by converging 
flows of incoming speed of 15 — 20 km s" 1 is fully compat- 
ible with Larson's relations. 

As no observational evidence for accretion onto 
T Tauri disks has been reported, our investigation of accre- 
tion driven turbulence in protostellar disks is more spec- 
ulative. However, very similar to galactic disks, without 
late mass accretion protostellar disks would drain onto 
their central stars on timcscales shorter than the inferred 
disk lifetimes. Using this as starting point, we were able to 
show that disk accretion can drive subsonic turbulence in 
T Tauri disks at roughly the right level if the rate at which 
gas falls onto the disk is comparable to the rate at which 
disk material accretes onto the central star. This process 
also provides a simple explanation for the observed rela- 
tion of accretion rate and stellar mass, M cx M*- 8 . 

We conclude that accretion-driven turbulence is a uni- 
versal concept with far-reaching implications for a wide 
range of astrophysical objects. 
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Appendix A: Motivation of the Potential Method 

The gravitational potential in rotationally supported disks is 
often approximately spherical symmetric. This holds for the ex- 
tended outer HI disks in dark matter dominated spiral galaxies 
(Section [3]) as well as for protostellar disks in the late phase of 
the evolution where most of the system mass is carried by the 
central star (Section [5]). The error is usually less than a few 
per cent. The equation of hydrostatic balance reads 



VP 



(A.l) 



with pressure P = c s p and gravitational force F g = 
GMp/ [R 2 + z 2 ) pointing towards the center. We use cylindri- 
cal coordinates, and c s and M are sound speed and enclosed 
mass, respectively. The pressure gradient in vertical direction 



dP 

dz 



H 



(A.2) 



where H is the vertical scale height. The z component of the 
force at any location H <C R is approximately 



F, 



GMp 
R 2 + H 2 

and consequently 

H c s 



H 



VR 2 + H< 



GMp H 
R 2 R 



R y/GM/R «rot ' 

with the circular velocity v Iot 



(A.l) 

(GM/R) 1/2 . This is Eq. Ip)l. 



Appendix B: The M- 
Densities 



cr 3 Relation at High 



In this appendix we discuss the validity of the relation eoc 
inferred in Section [2] 

First, we note from Fig. [2] that at high densities the veloc- 
ity dispersion of the dense gas remains nearly constant when 
the density threshold varies. This is clearly due to the fact 
that the total mass of the dense gas is constituted by several 
dense clumps which are not spatially correlated and randomly 
distributed in the turbulent box (as can be seen from Fig. [1]). 
Thus, the velocity dispersion simply reflects the velocity dis- 
persion of the box which typically varies with distance I as 
z i/3-i/2 xherefore] in this 

regime of density the integral of 
interest J pv 3 dV can be approximated as a 3 J pdV where a 
is approximately the velocity dispersion corresponding to the 
size of the more distant clumps. As a consequence, the depen- 
dence of Ma 3 on the density threshold should in this regime 
be identical to the one of the mass. 

For simplicity, let us consider isothermal gas. In 
this case, the density PDF is approximately log-normal 



o |= 
-l \ 

-3 \ 
-4 r 
-5 \ 
-6 L 




Fig. B.l. Fraction f t of mass above a threshold density 
Pt as function of density logarithmic density contrast St- 
For reference, the straight line indicates a slope of —1. 



(IVazauez-Semadenil Il994l : IPadoan et all 1 19971 : iKlessenl l2000l ; 
iKritsuk et al. mini), however, with higher-order correc- 
tions that become significa nt for highly compressive forcing 
(|Federrath et al.ll2008l . 120091 ). 



P{p) = 



■. exp 



(s-sf 



2a 2 



(B.l) 



where S = log(p/p), S = -a 2 /2 and cr 2 = ln(l + b 2 M 2 ) with 
parameter b ~ 0.5 and Alfvenic Mach number M ~ 5. 

The mass above some density threshold pt is obtained as 



Aft = pV 



%/27T<T" „ 



exp(<5) exp 



(S-Sf 



2a 2 



dS, 



where St = log (pt/p) and thus 



(A.3) f t = Mt/(pV) 



in I h 



exp (—a; 2 ) dx, 



(B.2) 



(B.3) 



-pV 



1 - erf 



S t - a 2 /2 
V2a 



Figure IBTT1 shows ft as function of St- Its behavior is sim- 
ilar to the mass distribution displayed in Fig. [2] and the re- 
lation Ma 3 oc p~ x is approximately recovered for St above 
~ 4 (where a is observed to be nearly constant). Thus, keep- 
ing in mind that the simulations include cooling, gravity and 
magnetic fields, we conclude that the effect we observe is rea- 
sonably well reproduced by this simple approach. At higher 
density contrast St, the slope of ft becomes steeper and we 
expect e to decrease more rapidly. 

References 

Agertz, O., Teyssier, R., & Moore, B. 2009, Monthly Notices 

RAS Letters, 397, L64, (c) Journal compilation © 2009 RAS 
Andre, P., Belloche, A., Motte, F., & Peretto, N. 2007, 

Astronomy and Astrophysics, 472, 519 
Andre, P., Ward- Thompson, D., & Barsony, M. 2000, 

Protostars and Planets IV (Book - Tucson: University of 

Arizona Press; eds Mannings, 59 
Audit, E. & Hennebelle, P. 2005, A&A, 433, 1 
Audit, E. & Hennebelle, P. 2010, Astronomy and Astrophysics, 

511, 76 

Avila-Reese, V. & Vazquez-Semadeni, E. 2001, ApJ, 553, 645 
Balbus, S. A. & Hawley, J. F. 1998, Reviews of Modern Physics, 
70, 1 



20 



Klessen & Hennebelle: Accretion-Driven Turbulence 



Ballesteros-Paredes, J., Gazol, A., Kim, J., et al. 2006, ApJ, 

637, 384, (c) 2006: The American Astronomical Society 
Ballesteros-Paredes, J., Hartmann, L., & Vazquez-Semadeni, 

E. 1999, ApJ, 527, 285 
Ballesteros-Paredes, J., Klessen, R. S., Low, M.-M. M., & 

Vazquez-Semadeni, E. 2007, Protostars and Planets V, 63 
Banerjee, R., Klessen, R. S., & Fendt, C. 2007, ApJ, 668, 1028, 

(c) 2007: The American Astronomical Society 
Banerjee, R., Vazquez-Semadeni, E., Hennebelle, P., & 

Klessen, R. S. 2009, Monthly Notices of the Royal 

Astronomical Society, 1130, (c) Journal compilation © 2009 

RAS 

Beck, R. 2007, A&A, 470, 539 

Beck, R., Brandenburg, A., Moss, D., Shukurov, A., & Sokoloff, 
D. 1996, Annu. Rev. Astro. Astrophys., 34, 155 

Bertout, C. & Genova, F. 2006, Astronomy and Astrophysics, 
460, 499 

Bigiel, P., Leroy, A., Walter, F., et al. 2008, The Astronomical 

Journal, 136, 2846 
Birnboim, Y. & Dekel, A. 2003, Monthly Notices of the Royal 

Astronomical Society, 345, 349, (c) 2003 RAS 
Blitz, L., Fukui, Y., Kawamura, A., et al. 2007, Protostars and 

Planets V, 81 

Blitz, L., Spergel, D. N., Teuben, P. J., Hartmann, D., & 
Burton, W. B. 1999, ApJ, 514, 818, (c) 1999: The American 
Astronomical Society 

Boldyrev, S. 2002, ApJ, 569, 841 

Bondi, H. 1952, Monthly Notices of the Royal Astronomical 

Society, 112, 195 
Bondi, H. & Hoyle, F. 1944, Monthly Notices of the Royal 

Astronomical Society, 104, 273 
Braun, R. & Thilker, D. A. 2004, A&A, 417, 421 
Brunt, C. M. 2003, ApJ, 583, 280 

Brunt, C. M., Heyer, M. H., & Low, M.-M. M. 2009, Astronomy 

and Astrophysics, 504, 883 
Calvet, N., Muzerolle, J., Briceno, C, et al. 2004, The 

Astronomical Journal, 128, 1294 
Casuso, E. & Beckman, J. E. 2004, A&A, 419, 181 
Ceverino, D., Dekel, A., & Bournaud, F. 2009, eprint arXiv, 

0907, 3271 

Chiappini, C, Renda, A., & Matteucci, F. 2002, A&A, 395, 
789 

Conselice, C. J. 2003, The Astrophysical Journal Supplement 
Series, 147, 1, (c) 2003: The American Astronomical Society 

Corbelli, E. & Salpeter, E. E. 1988, Astrophysical Journal, 326, 
551 

de Blok, W. J. G., Walter, F., Brinks, E., et al. 2008, The 

Astronomical Journal, 136, 2648 
Dekel, A. & Birnboim, Y. 2006, Monthly Notices of the Royal 

Astronomical Society, 368, 2 
Dekel, A., Birnboim, Y., Engel, G., et al. 2009, Nature, 457, 

451 

Dib, S. & Burkert, A. 2005, ApJ, 630, 238 

Dickey, J. M. & Lockman, F. J. 1990, IN: Annual review of 

astronomy and astrophysics. Vol. 28 (A91-28201 10-90). Palo 

Alto, 28, 215 

Dullemond, C. P., Natta, A., & Testi, L. 2006, ApJ, 645, L69 
Dutrey, A., Guilloteau, S., & Ho, P. 2007, Protostars and 
Planets V, 495 

Dziourkevitch, N., Elstner, D., & Riidiger, G. 2004, A&A, 423, 
L29 

Elmegreen, B. G. 2000, ApJ, 530, 277 

Elmegreen, B. G. & Burkert, A. 2009, eprint arXiv, 0912, 996 



Elmegreen, B. G., Elmegreen, D. M., Fernandez, M. X., & 

Lemonias, J. J. 2009a, ApJ, 692, 12 
Elmegreen, B. G. & Scalo, J. 2004, Annual Review of 

Astronomy &Astrophysics, 42, 211 
Elmegreen, D. M., Elmegreen, B. G., Marcus, M. T., et al. 

2009b, ApJ, 701, 306 
Elmegreen, D. M., Elmegreen, B. G., Rubin, D. S., & 

Schaffer, M. A. 2005, ApJ, 631, 85, (c) 2005: The American 

Astronomical Society 
Evans, N. J., Dunham, M. M., j0rgensen, J. K, et al. 2009, 

The Astrophysical Journal Supplement, 181, 321 
Falgarone, E., Hily-Blant, P., & Levrier, F. 2004, Astrophysics 

and Space Science, 292, 89 
Federrath, C, Duval, J., Klessen, R. S., Schmidt, W., & Low, 

M. M. M. 2009, eprint arXiv, 0910, 5469 
Federrath, C, Klessen, R. S., & Schmidt, W. 2008, ApJ, 688, 

L79, (c) 2008: The American Astronomical Society 
Ferriere, K. M. 2001, Reviews of Modern Physics, 73, 1031, (c) 

2001: The American Physical Society 
Ferrini, F., Marchesoni, F., & Vulpiani, A. 1983, Astrophysics 

and Space Science (ISSN 0004-640X), 96, 83 
Fich, M. & Tremaine, S. 1991, IN: Annual review of astronomy 

and astrophysics. Vol. 29 (A92-18081 05-90). Palo Alto, 29, 

409 

Field, G. B., Blackman, E. G., & Keto, E. R. 2008, Monthly 

Notices of the Royal Astronomical Society, 385, 181 
Fleck, R. C. 1983, Astrophysical Journal, 272, L45 
Fleck, R. C. 1996, Astrophysical Journal v.458, 458, 739 
Folini, D. & Walder, R. 2006, Astronomy and Astrophysics, 
459, 1 

Frisch, U. 1995, Turbulence. The legacy of A. N. Kolmogorov., 

price 45.00. ISBN 0-521-45713-0. Price 19.95. 
Frisch, U., Sulem, P.-L., & Nelkin, M. 1978, Journal of Fluid 

Mechanics, 87, 719 
Fromang, S., Hennebelle, P., & Teyssier, R. 2006, Astronomy 

and Astrophysics, 457, 371 
Fukui, Y., Kawamura, A., Wong, T., et al. 2009, ApJ, 705, 144 
Gatti, T., Natta, A., Randich, S., Testi, L., & Sacco, G. 2008, 

Astronomy and Astrophysics, 481, 423 
Gatti, T., Testi, L., Natta, A., Randich, S., & Muzerolle, J. 

2006, Astronomy and Astrophysics, 460, 547 
Genzel, R., Burkert, A., Bouche, N., et al. 2008, ApJ, 687, 59, 

(c) 2008: The American Astronomical Society 
Hartmann, L., Ballesteros-Paredes, J., & Bergin, E. A. 2001, 

ApJ, 562, 852 

Hartmann, L., D'Alessio, P., Calvet, N., & Muzerolle, J. 2006, 
ApJ, 648, 484 

Hartmann, L. & Kenyon, S. J. 1996, Annu. Rev. Astro. 

Astrophys., 34, 207 
Heiles, C. & Troland, T. H. 2005, ApJ, 624, 773, (c) 2005: The 

American Astronomical Society 
Heithausen, A., Bensch, F., Stutzki, J., Falgarone, E., & Panis, 

J. F. 1998, Astronomy and Astrophysics, 331, L65 
Heitsch, F., Burkert, A., Hartmann, L. W., Slyz, A. D., & 

Devriendt, J. E. G. 2005, ApJ, 633, L113 
Heitsch, F., Low, M.-M. M., & Klessen, R. S. 2001, ApJ, 547, 

280 

Heitsch, F. & Putman, M. E. 2009, ApJ, 698, 1485 

Heitsch, F., Slyz, A. D., Devriendt, J. E. G., & Burkert, A. 

2006a, Monthly Notices of the Royal Astronomical Society, 

373, 1379 

Heitsch, F., Slyz, A. D., Devriendt, J. E. G., Hartmann, L. W., 
& Burkert, A. 2006b, ApJ, 648, 1052 



Klessen & Hennebelle: Accretion-Driven Turbulence 



21 



Hennebelle, P. & Audit, E. 2007, Astronomy and Astrophysics, 
465, 431 

Hennebelle, P., Banerjee, R., Vazquez-Semadeni, E., Klessen, 

R. S., & Audit, E. 2008, A&A, 486, L43 
Hennebelle, P. & Chabrier, G. 2008, ApJ, 684, 395, (c) 2008: 

The American Astronomical Society 
Hopkins, A. M., McClure-Griffiths, N. M., & Gaensler, B. M. 

2008, ApJ, 682, L13 
Hunter, J. H., Sandford, M. T., Whitaker, R. W., & Klein, 

R. I. 1986, Astrophysical Journal, 305, 309 
Jones, B. F. & Herbig, G. H. 1979, Astronomical Journal, 84, 

1872, a&AA ID. AAA026.121.015 
Kalberla, P. M. W. 2003, ApJ, 588, 805, (c) 2003: The 

American Astronomical Society 
Kawamura, A., Mizuno, Y., Minamidani, T., et al. 2009, The 

Astrophysical Journal Supplement, 184, 1 
Klessen, R. S. 2000, ApJ, 535, 869 
Klessen, R. S. 2001, ApJ, 556, 837 

Klessen, R. S., Ballesteros-Paredes, J., Vazquez-Semadeni, E., 

& Duran-Rojas, C. 2005, ApJ, 620, 786 
Klessen, R. S., Heitsch, F., & Low, M.-M. M. 2000, ApJ, 535, 

887 

Kolmogorov, A. 1941, Dokl. Akad. Nauk SSSR, 30, 301 

Koyama, H. & Inutsuka, S. 2002, ApJ, 564, L97 

Kritsuk, A. G., Norman, M. L., Padoan, P., & Wagner, R. 

2007, ApJ, 665, 416 
Larson, R. B. 1981, Royal Astronomical Society, 194, 809, 

a&AA ID. AAA029. 131.011 
Larson, R. B. 2005, Monthly Notices RAS, 359, 211 
Leroy, A. K., Walter, P., Bigiel, P., et al. 2009, The 

Astronomical Journal, 137, 4670 
Leroy, A. K., Walter, F., Brinks, E., et al. 2008, The 

Astronomical Journal, 136, 2782 
Li, Y., Low, M.-M. M., & Klessen, R. S. 2005, ApJ, 620, L19, 

(c) 2005: The American Astronomical Society 
Li, Z.-Y. & Nakamura, F. 2006, ApJ, 640, L187 
Linsky, J. L. 2003, Space Sci Rev, 106, 49, (c) 2003: Kluwer 

Academic Publishers 
Lopez, R. G., Natta, A., Testi, L., & Habart, E. 2006, 

Astronomy and Astrophysics, 459, 837 
Low, M. M. 1999, ApJ, 524, 169 

Low, M. M. & Klessen, R. S. 2004, Reviews of Modern Physics, 

76, 125, (c) 2004: The American Physical Society 
Low, M. M., Klessen, R. S., Burkert, A., & Smith, M. D. 1998, 

Physical Review Letters, 80, 2754 
Lubowich, D. A., Pasachoff, J. M., Balonek, T. J., et al. 2000, 

Nature, 405, 1025, (c) 2000: Nature 
Madsen, S. 2003, Astronomy and Astrophysics, 401, 565 
Madsen, S., Dravins, D., & Lindegren, L. 2002, Astronomy and 

Astrophysics, 381, 446 
Makarov, V. V. 2006, The Astronomical Journal, 131, 2967 
Makarov, V. V. 2007, ApJ, 658, 480 
Matzner, C. D. 2002, ApJ, 566, 302 

Mckee, C. F. & Ostriker, E. C. 2007, Annual Review of 

Astronomy & Astrophysics, 45, 565 
Melioli, C, Brighenti, F., D'Ercole, A., & de Gouveia Dal Pino, 

E. M. 2008, Monthly Notices of the Royal Astronomical 

Society, 388, 573 
Melioli, C, Brighenti, F., D'Ercole, A., & de Gouveia Dal Pino, 

E. M. 2009, Monthly Notices of the Royal Astronomical 

Society, 399, 1089 
Minter, A. H. & Spangler, S. R. 1996, Astrophysical Journal 

v.458, 458, 194 



Mohanty, S., Jayawardhana, R., & Basri, G. 2005, ApJ, 626, 
498 

Muzerolle, J., Hillenbrand, L., Calvet, N., Briceno, C, & 

Hartmann, L. 2003, ApJ, 592, 266 
Muzerolle, J., Luhman, K. L., Briceno, C, Hartmann, L., & 

Calvet, N. 2005, ApJ, 625, 906 
Naab, T. & Ostriker, J. P. 2006, Monthly Notices RAS, 366, 

899 

Nakamura, F. & Li, Z.-Y. 2008, ApJ, 687, 354 

Natta, A., Testi, L., Muzerolle, J., et al. 2004, Astronomy and 

Astrophysics, 424, 603 
Natta, A., Testi, L., & Randich, S. 2006, Astronomy and 

Astrophysics, 452, 245 
Neuhaeuser, R., Sterzik, M. F., Schmitt, J. H. M. M., 

Wichmann, R., & Krautter, J. 1995, Astronomy and 

Astrophysics, 297, 391 
Ossenkopf, V. & Low, M.-M. M. 2002, A&A, 390, 307 
Ostriker, J. P. & Tinsley, B. M. 1975, Astrophysical Journal, 

201, L51, a&AA ID. AAA014.061.038 
Padoan, P., Kritsuk, A., Norman, M. L., & Nordlund, A. 2005, 

ApJ, 622, L61 
Padoan, P. & Nordlund, A. 1999, ApJ, 526, 279 
Padoan, P., Nordlund, A., & Jones, B. J. T. 1997, Monthly 

Notices of the Royal Astronomical Society, 288, 145 
Peek, J. E. G. 2009, ApJ, 698, 1429 

Piontek, R. A. & Ostriker, E. C. 2007, ApJ, 663, 183, (c) 2007: 

The American Astronomical Society 
Prochaska, J. X. & Wolfe, A. M. 2009, ApJ, 696, 1543 
Putman, M. E. 2006, ApJ, 645, 1164 
Roy, N. 2007, Monthly Notices RAS Letters, 378, L34 
Santillan, A., Sanchez-Salcedo, F. J., & Franco, J. 2007, ApJ, 

662, L19 

Scalo, J. & Elmegreen, B. G. 2004, Annual Review of 

Astronomy &Astrophysics, 42, 275 
Schmidt, W., Federrath, G, & Klessen, R. 2008, Physical 

Review Letters, 101, 194505, (c) 2008: The American 

Physical Society 
Sellwood, J. A. & Balbus, S. A. 1999, ApJ, 511, 660, (c) 1999: 

The American Astronomical Society 
Semelin, B. & Combes, F. 2005, A&A, 441, 55 
She, Z.-S. & Leveque, E. 1994, Physical Review Letters, 72, 

336 

Smith, M. C, Ruchti, G. R., Helmi, A., et al. 2007, Monthly 
Notices of the Royal Astronomical Society, 379, 755 

Spitoni, E., Matteucci, F., Recchi, S., Cescutti, G., & Pipino, 
A. 2009, Astronomy and Astrophysics, 504, 87 

Stone, J. M., Ostriker, E. C, & Gammie, C. F. 1998, ApJ, 508, 
L99 

Tamburro, D., Rix, H.-W., Leroy, A. K., et al. 2009, The 

Astronomical Journal, 137, 4424 
Teyssier, R. 2002, Astronomy and Astrophysics, 385, 337 
Thilker, D. A., Bianchi, L., Meurer, C, et al. 2007, The 

Astrophysical Journal Supplement Series, 173, 538, (c) 2007: 

The American Astronomical Society 
Throop, H. B. & Bally, J. 2008, The Astronomical Journal, 

135, 2380 

van Zee, L. & Bryant, J. 1999, The Astronomical Journal, 118, 
2172, (c) 1999: The American Astronomical Society 

Vazquez-Semadeni, E. 1994, Astrophysical Journal v. 423, 423, 
681 

Vazquez-Semadeni, E., Ballesteros-Paredes, J., & Klessen, 
R. S. 2003, ApJ, 585, L131 



22 



Klessen & Hennebelle: Accretion-Driven Turbulence 



Vazquez-Semadeni, E., Gomez, G. C, Jappsen, A. K., et al. 

2007, ApJ, 657, 870 
Vazquez-Semadeni, E., Gomez, G. C, Jappsen, A.-K., 

Ballesteros-Paredes, J., & Klessen, R. S. 2009, ApJ, 707, 

1023 

Vazquez-Semadeni, E., Ryu, D., Passot, T., Gonzalez, R. F., 

& Gazol, A. 2006, ApJ, 643, 245 
Vishniac, E. T. 1994, Astrophysical Journal, 428, 186 
von Weizsacker, C. F. 1951, Astrophysical Journal, 114, 165 
Wakker, B. P., Howk, J. C, Savage, B. D., et al. 1999, Nature, 

402, 388, (c) 1999: Nature 
Walder, R. & Folini, D. 1998, Astronomy and Astrophysics, 

330, L21 

Walder, R. & Folini, D. 2000, Astrophysics and Space Science, 
274, 343 

Walter, F., Brinks, E., de Blok, W. J. G., et al. 2008, The 

Astronomical Journal, 136, 2563 
Wang, P., Li, Z.-Y., Abel, T., & Nakamura, F. 2010, ApJ, 709, 

27 

Wichmann, R., Krautter, J., Schmitt, J. H. M. M., et al. 1996, 

VizieR On-line Data Catalog, 331, 20439 
Xue, X. X., Rix, H. W., Zhao, G., et al. 2008, ApJ, 684, 1143, 

(c) 2008: The American Astronomical Society 
Zaritsky, D. & Christlein, D. 2007, The Astronomical Journal, 

134, 135 



